arXiv:1509.05695v3 [math.AP] 27 Apr 2016 


SOLVABILITY OF THE DIRICHLET, NEUMANN AND THE REGULARITY 
PROBLEMS FOR PARABOLIC EQUATIONS WITH HOLDER CONTINUOUS 

COEFFICIENTS 

ALEJANDRO J. CASTRO, SALVADOR RODRIGUEZ-LOPEZ, AND WOLFGANG STAUBACH 


Abstract. We establish the L^-solvability of Dirichlet, Neumann and regularity problems for 
divergence-form heat (or diffusion) equations with time-independent Holder-continuous diffusion 
coefficients, on bounded Lipschitz domains in R”. This is achieved through the demonstration 
of invertibility of the relevant layer-potentials which is in turn based on Fredholm theory and a 
systematic transference scheme which yields suitable parabolic Rellich-type estimates. 
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1. Introduction 

In this paper, we prove the solvability of the Dirichlet, Neumann and Regularity problems (DNR 
problems for short) for divergence-form parabolic equations of the form 

dMX,t) - Vx • (A(A)Vxu(V,t)) = 0 

on bounded Lipschitz domains in K" (n > 3), under the assumptions that A{X) is uniformly elliptic, 
symmetric and Holder-continuous. 
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Let US very briefly recall some of the basic results in the held of second-order parabolic boundary 
value problems with time-independent diffusion coefficients on low regularity domains, which are 
the predecessors of the present paper. For the sake of brevity, we confine ourselves to only mention 
those investigations for parabolic equations which have dealt with the solvability of the boundary 
value problems mentioned above. In [^, E. Fabes and N. Riviere proved the solvability of the 
Dirichlet and Neumann problems on bounded domains. This paper paved the way for subsequent 
developments in the held. Since the appearance of [S], the investigations of the parabolic boundary 
value problems have been concerned with either lowering the regularity of the boundary of the 
domain, or lowering the regularity of the matrix A appearing in the equation, or both. Another 
goal is to consider U’ boundary value problems for various values of p (i.e. p other than 2). 

But in this paper, it is the boundary value problems for time-independent domains that concern 
us. In [5] E. Fabes and S. Salsa investigated the caloric measure and the {p > 2) solvability 
of the initial-Dirichlet problem for the usual heat equation in Lipschitz cylinders. In paper [2], R. 
Brown studied the boundary value problems and the layer potentials for the heat equation on 
bounded Lipschitz domains in K". Next step was taken by M. Mitrea in m where he proved the 
LP solvability (for suitable values of p), of DNR problems for divergence-type heat equations with 
smooth diffusion coefficients on compact manifolds with Lipschitz boundary. Our investigation in 
this paper is the continuation of these lines of studies by further pushing down the regularity of the 
diffusion coefficients and assuming only Holder continuity, which together with the assumptions of 
ellipticity and symmetry, will yield the solvability of DNR problems. 

In m, M. Mitrea and M. Taylor proved the solvability of the DNR problems for elliptic equations 
involving the Laplace-Beltarmi operator with Holder-continuous metrics on Riemannian manifolds 
with Lipschitz boundary. Later in [S], C. Kenig and Z. Shen used the method of layer potentials to 
study boundary value problems in a bounded Lipschitz domain in K", with n > 3, for a family 
of second-order elliptic systems with rapidly oscillating periodic coefficients. As a consequence, 
they also established the solvability of the DNR problems for divergence-form elliptic equations 
Vx • {A{X)Wxu{X)) = 0 on the aforementioned domains, under the assumptions that A{X) is 
uniformly elliptic, symmetric, periodic and Holder-continuous. Our paper could be considered as a 
parabolic counterpart of [8] and HU. 

We shall now briefly describe the main results of the paper and the structure of this manuscript. 
We recall that a bounded domain D C K” is called a Lipschitz domain (with Lipschitz constant 
M > 0) if dLt can be covered by finitely many open circular cylinders whose bases have positive 
distance from dil, and corresponding to each cylinder Z C K" there exists: 

• a coordinate system (x',Xn), with x' G and G R, such that the a;„-axis is parallel 

to the axis of Z\ 

• a function p : R"”^ —R satisfying the Lipschitz condition 

\p{x')-p{v')\<M\x' -y'\, x',y'GR”-\ 

such that 

(1) LtC\ Z = {{x\ Xn) G Z : Xn > ‘p(x')} and dfi Cl Z = {(a;', a;„) G Z : Xn = p{x')}. 
As mentioned earlier, we consider the parabolic divergence-type equation 

(2) Cau{X, t) := dtu{X, t)-Vx- (A(A)Vxw(A, t)) = 0 in x (0, T), 

where 0 < T < oo and LI is an open bounded Lipschitz domain in R", n > 3. We assume that the 
real matrix A{X) = {aij{X)) verifies the following properties: 

(Al) Independenee of the time-variable: A = A(A); 

(A2) Symmetry: Gij = ap, i, j = 1,... ,n; 

(A3) Uniform ellipticity: for certain /i > 0, 

n 

MieP < E < -ICP, G R’^; 

(A4) Holder regularity: for some k > 0 and 0 < a < 1, 

\a,jiX)-a,j{Y)\<K\X-Yf, X,YgR^, i,j = l,...,n. 
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To state the aforementioned DNR problems, one defines the lateral boundary of D x (0, T) as 
St ■= dn X (0,T). Moreover the conormal derivative associated with the operator Ca will be 
defined as 

(3) d^u{Q,t) := d^j,u{Q,t) := {VYu{Y,t)\^^Q,A[Q)NQ), {Q,t) & St, 

where Nq = (ni((5),... ,nn{Q)) denotes the unit inner normal to dfl at Q, which is defined a.e. 
on dfl. The conormal derivative is sometimes denoted by 9,^^ to emphasise its dependence on the 
matrix A. One can also define the tangential derivative Vy of a function u by 

( 4 ) VTu(Q,t) := - (Vvu(y,t)|,^Q, 7 VQ)iVQ, {Q,t) G St- 


Given these preliminaries, we are interested in the solvability, in the weak sense, (see Section]^ for 
the proper statements) of the following problems: 


Dirichlet’s problem 


Neumann’s problem 


{ Cau = 0 in f2x(0,r) 

u(X,0) = 0, XGfl 

u = f G L‘^{St) on St 


{ Cau = 0 in Ox (0, T) 

u{X,0)=0, XGn 
dyU = f G L'^{St) on St 


Regularity problem 

{ Cau = 0 in Ox (0,T) 

u(X,0) = 0, XGfl 
u = / e on St 


To achieve our goals, in Section we introduce the notations and recall some basic harmonic an¬ 
alytic tools which will be used throughout this paper. In Section we prove quite a few new 
estimates for various derivatives of the fundamental solution of parabolic divergence-type opera¬ 
tors with Holder-continuous diffusion coefficients, and also prove the corresponding estimates for 
the Fourier transform of the fundamental solution in the time variable. It should be noted that 
although the estimates that are obtained are similar to those in the constant coefficient case, this 
doesn’t simplify the study of the solvability of the DNR problems in our setting. Indeed even in 
the elliptic divergence-form case studied in |8] and m, one also has the same estimates as those 
for the constant coefficient Laplacian, but that by no means simplifies the problem. The major 
difficulty in the study of low regularity elliptic and parabolic problems is to show the invertibility 
of the corresponding layer potentials which is a significant task for equations with rough coefficients. 

In Section]^ we study the parabolic single and double-layer potentials associated to the operator Ca 
and establish the boundedness of the boundary singular integral corresponding to this operator 
as well as the boundedness of the non-tangential maximal function associated to the double-layer 
potential and that of the normal derivative of the single layer potential. In this section we also prove 
a couple of jump formulas for the aforementioned operators and also show the boundedness of the 
non-tangential maximal function associated to a fractional derivative of the single layer potential. 
The proofs of the boundedness here is somewhat simpler since it can be done using parabolic 
Calderon-Zygmund theory, as was carried out by Brown [5] in the constant coefficient case. Next, in 
Section]^ we consider the Fourier transform of the equation CAu{X,t) = 0 (in time) and establish 
the estimates proved in Section for the Fourier-transformed equation. These estimates will be 
very useful for us in one of the central sections of our paper namely Section 

Here we have an approach which allows us, to transfer in a systematic way, estimates for equa¬ 
tions with smooth coefficients to those with Holder-continuous coefficients. Briefly, the transference 
method works as follows. One writes the original Holder-continuous diffusion matrix A as B + C 
where B = A + A — and C = — A, with a smooth diffusion coefficient A and a suitable A^’’^ 
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such that IICIloo = — 4||oo can be made arbitrary small by choosing r small enough. Then 

the first step is to prove Rellich estimates for the smooth part A and then transfer those estimates 
to B. Moreover, those terms in the invertibility estimates for the operator associated with A that 
involve C, can be handled using the smallness of ||C||oo and suitable boundedness estimates. 
Apart from the proofs of the solvability of DNR problems, this transference method is one of 
the main achievements of the present paper. 

In Section [7] we prove the invertibility results which are the key to the solvability of the DNR 
problems. This is done by using all the information that we have gathered up to that point and 
an application of Fredholm theory. Finally in Section we very briefly outline the solvability of 
the Dirichlet, Neumann and regularity problems, which is as usual, a standard consequence of the 
invertibility of the relevant singular integral operators. 

Acknowledgements. The authors would like to thank Kaj Nystrom for bringing the subject of low 
regularity parabolic boundary value problems to their attention and for stimulating and encouraging 
discussions on this topic. We are also grateful to Carlos Kenig for his encouragement and interest 
in our work. 


2. Basic notations and tools 

One of the conventions in the theory of boundary value problems on low-regularity domains, which 
we shall follow hereafter, is that interior points in the domain will be denoted by A, Y while those 
of do will be denoted by P,Q. Furthermore, it is also important to warn the reader that, when we 
write dP or dQ in the integrals that are performed over the boundary, then dP or dQ denote the 
surface measures da{P) or da{Q). 

We sometimes write a < 6 as shorthand notation for a < Cb. The constant C hidden in the estimate 
a Y b can be determined by known parameters in a given situation, but in general the values of such 
constants are not crucial to the problem at hand. Moreover the value of C may differ from line to line. 

Now let D’*' := and n~ := M”\D. Then for some a > 0 one defines the non-tangential approaching 
domains 7+(F’) C D+ and 7 -(P) C as follows: 

(5) 7 +(P) := {X €n:\X - P\ < (1 -h a) dist(A, 911)}, 

(6) 7 _(P) := {A G K”\n: |A-P| < (1-h a) dist(A,. 

It is important to note that, for every P,Q G dil and A G 7±(P) one has 

(7) |A-Q| >dist(A,9D) > :^|A-P|, 

1 -|- a 

and 

|A - Q| > |P - Q| - |A - P| > |P - Q| - (1 + a)|A - Q\ 

which implies, 

(8) |X-Q|>^|P-Q|. 

2 -I- a 

These estimates will be used in Sections]^ andin connection to the estimates for non-tangential 
maximal functions associated to various operators. 

Given ^ and (|^ and a function u, for eve 
is defined by 

(9) ut{P,t) 

We consider also the non-tangential limits 

(10) u+iP,t) 

-(PA) 


3ry (P, t) G St, the non-tangential maximalfunction u* 


:= sup \u{X,t)\- 


:= lim |tt(A, f)|, 
x^p 
xej+{P) 


:= lim |it(A, t)|. 

X-s-P 

xe-/-(P) 


( 11 ) 


u 
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The two limits defined above are the ones that appear in the jump relations occurring in this paper, 
see Sections Hand [3 


We denote by and Aion the Hardy-Littlewood maximal operators on R and dQ respectively, 
that is, 

Mi{h){t) = sup ^ [ |h(s)|ds, 

r>0 J|t-s|<r 


and 


MaQ{g){P) = sup^^ 

r>0 T 


I aQn{\P-Q\<r} 


\g{Q)\dQ. 


It is well-known that both operators are bounded in L^. We write M to refer to A4i or Mon, 
indistinctly. We also recall the following well-known result which will be used in the proofs of our 
estimates in Sections 0 and H 


Lemma 2.1. Let (j) be a positive, radial, decreasing and integrable function. Then 

sup \(j3r * < \miMiF){u}), 

r>0 

where ^ is the usual convolution in fr) jr'^ and m = 1 or m = n — 1. 

We shall also make a repeated use of the so called Schur’s lemma 


Lemma 2.2. Let X, Y be two measurable spaces. Let T be an integral operator with Schwartz 
kernel K{x,y), x & X, y &Y 


If 


for almost all x and 


Tfix) = J K{x,y)f{y)dy. 

J \K{x,y)\dy<a 

[ \K{x,y)\dx < 13 


JX 


for almost all y, then T extends to a hounded operator T : L‘^{Y) L‘^{X) with the operator norm 


\\T\\l^{y)^l^{x) < y/chS- 

For the application of the Fredholm theory in Section H we would also need the following elementary 
functional analytic lemmas. We include the proofs for the convenience of the reader. 


Lemma 2.3. Let 6>0, n be a bounded domain in R" and C be the operator defined by 

Then, C is a compact operator in L‘^[dLl). 

Proof. We write 

C{g){P) = [ K{Q- P)g{Q)dQ, P e 90, 

Jon 

where K{Z) := Analogously we consider, for each e > 0, the operator Cg associated to 

the kernel Kg{Z) := {\Z\ Since, 

for each e > 0, Cg is a Hilbert-Schmidt operator and hence it is compact. Moreover Lemma [2T^ and 
the Lebesgue dominated convergence theorem yield 

ne'e — C\\L'^(dn)^L^{dn) < \\Ke — Ar|jLi(an) —0, as e —>• 0. 

Therefore C is a compact operator. □ 


In the estimates for the difference of the parabolic single layer potentials associated to two different 
diffusion coefficients, the following equality from the theory of Markov chains is useful. 
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Proposition 2.4. {Chapman-Kolmogorov formula) Let r, s and t be real numbers with r < s <t 
and let A > 0. Then 


( 12 ) 


/ 






g 2(t-s) g 2 (,-t-) 




1'^/t — r\~'^/'^ -A|™--u| 


g 2(t-r-) 


See 13 Proposition 3.2.3] for a proof of the Chapman-Kolmogorov formula. 


We shall also need the following elementary functional analytic lemma, which is useful in connection 
to the invertibility of the boundary singular integral: 

Lemma 2.5. Let T be a bounded linear operator from a Hilbert space H into itself. Furthermore, 
assume that T is injective, has closed range, and that T* — T is a compact operator. Then T is also 
surjective. 

Proof. Since T is a bounded, injective and has closed range, it is well-known that ind (T -|- K) = 
ind (T) for all compact operators K : H ^ H, where indT := dim Ker T — dim Coker T, denotes the 
Fredholm index of T. Therefore, since ind(r) = —ind(T*) and ind {T*) = ind {T-\-T* —T) = ind (T), 
due to the compactness assumption on T* — T, we also have that ind (T) = —ind(T). Therefore 
ind (T) = 0, which together with the injectivity of T, yields the surjectivity of T. 

a 


Lemmas 2.6 and 2.7 below are taken from jS] and will be used in the proof of the Parabolic Rellich 
estimates in Section |6l 

Let ro := diam Ll < oo. We choose a cube Qa C K” such that ft c Qn- We call 2Qq to the cube 
with the same centre as Qq but with the double size. 


Lemma 2.6. [HI Lemma 7.1] Given a matrix A sa tisfyin g prop erties 1(A1)| - 1(A4)[ there exists 
A S C°°(2Qn \ dLl), such that A = A on dLl. Moreover (41) - (44) hold for A with a certain Holder 
exponent ao G (0, a] and 

(13) s diS(W^' 


Lemma 2.7. jS] Lemma 7.2] Let A be a matrix satisfying properties [(41)] - [(44)] Fix 9 G 


C')?°(—2ro) 2ro) such that 0 <9 <1 and 9 = 1 on (—ro,ro). Define, for each 0 < r < 1, 


4^(4:) := 


3dist(X, dfl) 


A{X), X G 2Qn, 

with the same 


where A is the matrix given in Lemma |2.6| Th en 4*^’’^ satisfies properties (41) - (44) 
Holder exponent oq € (0,a] as in Lemma 2.6 Moreover ||4(’') — 4||oo ^ 


In the investigation of the solvability of the regularity problem, we would need to deal with fractional 
Sobolev spaces. The fractional derivatives are defined as follows: 

Let / G C°°(—oo,T) and f(t) = 0 for t < 0. Then letting r(cr) denote Euler’s gamma function, one 
defines the fractional integrals and fractional derivatives of / via 

and 


(14) 


Dffit) 


dJi-afit), for 0 < cr < 1 
dtf{t), for cr = 1. 


Furthermore, for cti, CT 2 in (0,1) and cti -|-172 < 1 one has the following identities for the fractional 
integrals and derivatives 


DA{DTU))=DA^'^^f. 


{2ttY 

V2 


(1 -kisign(T)) |T|'"/(r), 


Note that one also has 
(15) 


Dffir) 
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where /(r) = dt. This will be useful in connection to the estimates for the fractional 

derivative of the single layer potential. 

Now given St as in the introduction section of this paper, the fractional Sobolev space H^’^^‘^{St) is 
the closure of space {n; v = uls^,, u G u{X,t) = 0 for t < 0} with respect to the norm 

■■={ r [ + Inp) dPdt 

Jo Jdn 



1/2 • ... . I I 

where ' is the fractional derivative defined using (141 and Vt is the tangential derivative defined 

in Q. 

The following estimate involving fractional derivatives, which has been taken from Brown’s thesis 


|5], will play an important role in the proof of parabolic Rellich estimates in Subsection 6.2 


Lemma 2.8. Let f,g G C°°{—oo,T) and f{t) = g{t) = 0 for t < 0. Then, 

1 / 2 , 


a i \ 1/2 / 

\f{t)?dtj \Dy*{g){t)\^dt 


1/2 


We conclude this section by pointing out that in what follows, due to the elementary and standard 
nature of the arguments and lack of space, we will follow the common practice of refraining from 
comments on justifications of the legitimacy of interchanging the order of integrations and that of 
differentiations and integrations. Certainly, all these operations can be fully justified in each case 
under consideration by a careful glance at the relevant proofs. 


3. Estimates for the fundamental solution of Ca with a Holder-continuous matrix 
Let r and T* be the fundamental solutions in K" for the operators Ca and C\ respectively, that is 
CaT{X, t- Y, s) = S(X - Y)S(t - s), C*aT*{Y, s; X, t) = 5{X - Y)5{t - s), 

where 5 denotes the Dirac’s delta function and C\ is the adjoint of Ca- Note that, if A is symmetric, 
C\u = —dt — V • (HVm). Also one has 

( 16 ) r(A,t;y,s) = r*(y,s;A,t), X,YGW^,t^s, 

see e.g. |3l Lemma 3.5]. Moreover, in the case of time-independent matrix A, we have that 
r(A,t;r,s) = r(A,t-s;r,o), t^s. 

From now on, we simply write the three-argument function T{X,Y,t — s) to refer to the above 
quantity, when there is no cause for confusion. 

Recall the relation 

r{x,Y,t)dt = f{x,Y), A,yeK”, 

where T represents the fundamental solution to the elliptic operator V • (AV-), see e.g. [TJ p. 895]. 
Next we collect some pointwise estimates for the fundamental solution, which shall play a basic role 
for the estimates of the forthcoming sections. Note that the constants appearing in the estimates 
below will depend on various combinations of n, fi, k and a, see the introduction section for the 
definitions of these latter constants. 



The following two lemmas are well-known for the fundamental solutions of divergence-type operators 
under much weaker conditions than those stated here, but since the regularities lower than Holder- 
continuity don’t concern us in this paper, we confine ourselves to the statements below. 


Lemma 3.1 (|T]). Assume f/taf ](A3)] holds. 

\T{X,f,Y,s)\< 


Then for every X,Y G K" and t, s > 0 one has that 

g-c\X-Yf/(t-s) 

_ s)n/2 (^)- 


Lemma 3.2 (' |12l Property 10, p. 163]). Assume that \{Al) and (A3) hold. Then, for every £ £ N, 
X,Y G K" and t > 0 one has 

-c\X-Y\^/t 

l^tr(A,y, t)| < X(o,oo)(/)- 


Later on, in proving the estimates we would also need estimates on the spatial derivatives of the 
fundamental solution. 
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Lemma 3.3. Assume that 
and t > 0 we have 


(Al) 


(A4) hold. 


Then for every m G N" such that \m\ < 2 ; X,Y G M" 


\d^r{x,Y,t)\ + \dipnx,Y,t)\ < 

where .9™" if m = (mi,..., m„) G N" and |m| = mi + • • • + m„. 

Proof. The bound for X-derivatives can be found in [HI eq. (13.1), p. 376]. To obtain the corre¬ 
sponding estimate for the Y derivative, define r{Y,X,t) := T* (Y, X, —t) = r{X,Y,t) and observe 
that it satisfies £yir(F, X, t) = 6{X — Y)S{t). Therefore 

\d^r{x,Y,t)\ = \dtpr*{Y,x,-t)\ = \d^TiY,x,t)\ < ^-^^l^x(o,oo)(i). 

□ 


The following few lemmas are entirely new and will be useful in the later sections. 


Lemma 3.4. Assume that (Tl)| - (T4)| hold, 
that 


Then, for every £ G N; X, y G K” and t > 0 we have 


c|X-Y|Vt 

(i) \dxdYTiX,Y,t)\ < ■■■ ■■X{0,oo)it), 

ill) \dxdiV{X,Y,t)\ + \dYdlT{X,Y,t)\ < 

where dx = dx^, for some j = 1,... ,n. 


Proof. We only prove (i) in detail and briefly comment on the proof of (ii). Fix X,Y G M", 
0 < s < t. Using the well-known semigroup property of the fundamental solution (see e.g. (2.41) in 
0), namely 

T{X,Y,t-s)= I r{X,Z,t-r)r{Z,Y,r - s)dZ, 


for any r G (s, t), we have 


dxdYT{X,Y,t- s) = [ dxT{X,Z,t-r)dYT{Z,Y,r- s)dZ. 

Thus the estimates in Lemma |3.3| and Proposition |2.4| yield 

\dxdYriX,Y,t- s)| < ((t - r)(r - s)(t - s)")-i/2e-d^-^lVb-s). 

Taking r = {t + s)/2 we obtain (i). 

To prove {ii), we observe that 9t(r(X, Y,t — s)) = —i9s(r(X, Y,t — s)), so following the same line of 
argument as in (z) using lemmas 3.2 and 3.4 yield {ii). 

□ 


For the fractional derivative defined as (141, we have the following estimates: 

Lemma 3.5. Assume that (Al) - (A4) hold. Then, for every X,Y G K" and t > 0 we have that 

(*) \ Dy ^ nx , Y , t )\ < ^3;,|^''_y|l, x(o.oo)(£), 

(zz) \dxDl/^T{X,Y,t)\ + \dYDl^^r{X,Y,t)\ < (t), 

(zzz) \ dtDy ^ r { X , Y , t )\ < -^^^^^y(0.oo)(£). 


Proof. To prove (z) it is enough to assume that X Y, which according to Lemma 3.4 part {ii) 
yields the continuity of r(X, Y,t) for X ^ Y and t G ( 0 , 00 ). Therefore using the continuity of 
T{X,Y,t), the definition of the fractional derivative (with r(l/ 2 ) = y/n), and integration by parts, 
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we obtain 


y^Dy^T{X,Y,t) = dt [ 
Jo 


T{X,Y,t-s) ^_ T{X,Y,0) ^ dtT{X,Y,t-s) 


r(x,r, 0 ) 

Vi 


+ 


Vi 

dsT{X,Y,s) 


ds 


— . 


ds 


T{x,Y,t/2) r(x,y,s) 




lo 2(t-s)3/2 


ds 


I* 9.r(x,y,s) 
It/2 y/t — s 
/■* a,r(x,y,s) 


/t/2 


Vt - 


ds 

' ds. 


Now using Lemma 3.2 to estimate each of the three terms above yields (i). The proofs of {ii) and 
(in) differ marginally from that of (i), however in proof of (ii), instead of using Lemma |3.2| one 
has to use Lemma 3.4 {ii). The details are left to the reader. □ 

In our transference scheme which would enable us to transfer invertibility of layer potential operators 
associated to smooth diffusion coefficients to the invertibility of non-smooth layer potentials, the 
following simple lemma is very useful. 

Lemma 3.6. Let Ai and A 2 be two diffusion coefficients, with the corresponding fundamental 
solutions (X, Y,t — s) and T^^ Y,t — s). Then the following equality holds for the difference 
of fundamental solutions: 

rAffx,Y,t-s)-rAVx,Y,t-s) 

^ pOC) p 

= V/ / d,YAAX,Z,u)d,YAVZ,Y,t-s-u){AffZ)-A2{Z))dZdu. 

^,J=l -^0 "'R" 

Proof. Integration by parts yields 

^ pOO p 

Y. / / d,TAffX,Z,u)d,YAVZ,Y,t-s-u)iAi{Z)-A2{Z))dZdu 

-^0 

Tl poo p 

= Yi / d,ffA 2 {Z)d,TAVZ,Y,t-s-u))rAAX,Z,u) dZdu 

^ Jo Jr- 

- [ [ dzV^iiZ)dziPAAZ,Y,u))TAiiZ,Y,t - s - u) dZ duy 

Now the claimed equality follows by also observing that 

dtVA^ {Z,Y,t- s-u)-X ■ {A 2 {Z)XV aVZ, Y,t-s-u))= 5{Y - Z)5{t-s- u), 
duTA^ {X, Z, u) - V • (All iZ)XTA, (X, Z,u))=S(X-Z) S(u), 
and that 

du [ TAAX,Z,u)dZ = dt ( TaAZ^YA- s-u)dZ = Q, 

Jw- Jr- 

since the integrals that are being differentiated are both equal to 1 regardless of the time variable. □ 

In [S] , the problem of the invertibility of boundary singular integrals was handled by utilising the 
time independence of the Laplacian in the heat equation and performing a Fourier transformation 
in the time variable. This is an approach which we also adapt here and it has numerous advantages. 
However, it behoves us then to get suitable estimates for the fundamental solution of the Fourier- 
transformed operator. To this end, we define the truncated Fourier transform of a function h as 


hij) ■.= Xt{h){T) := 


*'h{t) dt. 


If {Al) is satisfied, we can take the Fourier transform in time in and get the new equation 

(17) Cau{X, t) := -iTu{X, r) - Vx • {A{X)Vxu{X, r)) =0, X e fl, 

for each t. This way, the parabolic equation becomes an elliptic equation depending on the param¬ 
eter T, which we assume to be fixed hereafter. Moreover, it is clear that 

poo 

f(X,y,T)= / e-^^V{X,Y,t)dt, X,yGK’", 


is the fundamental solution of The following lemmas establish estimates for r(X, F, r). 
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Lemma 3.7. Assume that (Al) and [(.43)] hold. Then, for every N gN; X,Y G M" we have that 

|r(.Y,r,T)|<- \x-Y\^-^ -’ 

Proof. An integration by parts and Lemma |3.2| lead to 

Y, \x - yps) ds 


i\T\\iX -Y)\T\f{X,Y,T)\ = \X -Y\^\ 9f(, 

nOO nOC) 

<\X-Yf \dfr{x,Y,\x-Yfs)\ds = \x-Yf+^^ |(9f r)(A, y, |x - 

Jo Jo 

< |A-y|2+2^ f 
Jo 


ds 


=-c/s 


ds 5, 


(|A-y|2s)"/2+^ ^|A-y|’"-2' 

For various derivatives of F one also has the following estimates: 


□ 


Lemma 3.8. Assume that (Al) - (A4) hold. Then for every q > 0, m G N” such that \m\ < 2 and 
X,Y G K" we have that 

(*) |9^f(A,y,r)| + |5??r(A,y,T)|< ^ 


{a) |i9x9Yr(A,y, r)| 


< 


1 


|A-y|"’ 

(m) |9xr(A,y,r)| + |avf(A,y,r)|< 


{iv) |avr(A,y,n) - dyVix, y, t 2 )| 


< 


|A-y|"-2+" 


|A-y|"-i+29 

In - t 2 |^ 
|A-y|"-i-2/5 


, 'T 0, 

, for n > 3 and all fi G (0,1). 


Proof, (i) and (ii) are straightforward applications of Lemmas 3.3 
For (in), if g = iV G N, we can proceed as in the proof of Lemma 


and 3.4 (i). 

3.7 taking into account Lemma 


|3.4| (a). Finally, if q is not an integer, we use the following simple interpolation argument. Namely, 
write q = N + 9, where fV = [gj G N and 6* G (0,1). Then using (in) for the integer values of g, we 
obtain 

|axr(A,y,r)| = \dxfix,Y,T)\^-^\dxfix,Y,T)f 
< 


f k|-^ ^ 

1-0/ |t-|-(JV+1) 

I® |r|-« 

v|A-y|"-i+2^/ 

) V|x - y|"-i+2(A®+i), 

1 Jj^_yn-l +29 


The proof of the estimate for SyT is exactly the same. 

Statement (iv) is a consequence of the elementary estimate | < |f(ri — T 2 )|^, valid for 

all 0 < /? < 1 and Lemma [3. 3 1 Indeed we have 


|9yr(A,y,n) - dyVix, y, t2)| < 


—itri —itT2 


-e-^^^^\\dyr{X,Y,t)\dt 


poo poo 

<|Ti-'r 2 |'^/ t'^|i9yr(A,y,t)| dt < In - T 2 |^ / t^^ 

Jo Jo 


oo ^g_c|X-Y|Vt 


f{n+l)/2 


— dt 


< 


ri 


- T2|^ 


|x-y|"-i-2/3 


g-s^(n-3-2/3)/2 < 




- T2\^ 


\X - y|n-l-2/3 ’ 


□ 


provided that n > 3 and (3 G (0,1). 

For the Rellich estimates in Section we would also need the following general lemma: 

Lemma 3.9. Assume that (Al) - |(A4)| hold. Let g = 0 or g = 1/2 and let B be the operator defined 
by 

B{g){P):= f \T\^f{P,Q,T)g{Q)dQ, P G dLl. 

Jon 

Then, 

\\Bg\\L'^(dQ.) ^ IlfflliRon): 9 € 
where the estimate is uniform in r. 
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Proof. By Lemma [2. 2 1 and the symmetry of the kernel, it is enough to show that 



|r«r(P,Q,T)|dQ<l, 


PGdn, 


uniformly in r. 

If q = 0, by Lemma 3.8 (i), we just need to check that 


f dQ 

Jan IQ - P\^-^ 


p e dil. 


Locally, we can write P = {P',ip{P')), Q = {Q'P',Q' G for a certain Lipschitz 

function (p. Moreover, since ||Vvj||loo < 1 and is a compact set, there exists M > 0 such that 
IQ' — P'\ < M, for every Q,P G did. Then, the above integral is equal to 




Q'GK"-i (IQ' - P'|2 + - (p(p/)|2)(n-2)/2 


\Q'-P'\<M 


dQ' < /,. 


dQ' 


gb"-i |n'_p'|"-2 

\Q'-P'\<M 


< 1, P G dVL. 


Suppose now that q = 1/2. We are going to proceed as before but with a slight modification in 
order to avoid the dependence on the parameter r. We use the improved estimate in Lemma |3.7| 
with iV > 1, and write the corresponding integral in as 



- (IQ, + -VI + IVpW )F dQ 


< 


f |r|V^min{l,(|r||Q'-P'p)-^} 

/q-GB"- |Q'-P'h-2 


min{l, |Z|-2^} 


--in{|z|<i} 


dZ 

Zn-2 


dZ 


, ,, , , 7n-2+2Af 

■*-in{|z|>i} ^ 


VZGB’*-i 
< 1, P G dil. 


dZ 


□ 


The following lemma will also be useful in dealing with the transference of the invertibility of 
boundary singular integrals. 


Lemma 3.10. Let Ai and A 2 be two diffusion coefficients, with the corresponding fundamental 
solutions TAi{X,Y,t) and Ta 2 {X,Y,t). Then the following equalit'y holds for the difference of 
fundamental solutions: 



a,, f A, {X, z, t) Ta, {Z, Y, t) (VI2 (Z) 


Ai{Z))dZ. 


Proof. The proof is a consequence of Lemma |3.6| and taking the Fourier transform in the time 
variable. □ 


4. Parabolic layer potential operators; SLP, DLP and BSI 


The main operators, concerning elliptic and parabolic boundary value problems are the Layer- 
potential operators. One defines the parabolic single and double-layer potential operators by 


and 


S{f){X,t) 


nx,Q,t 


0 JdQ. 


s)f{Q,s)dQds, X Gil, t > 0, 


'Pif)iX,t) := [ [ d,2TiX,Q,t-s)f{Q,s)dQds 
Jo Jan 

pt p 

= X! / / a^ffQ)nj{Q)dyT{X,Y,t- 

A A — t Jo JdO. 




=-.Y,V^'={f){X,t), XGil,t>Q. 

* j=i 


s)w^Qf[Q,s)dQ ds 


The single and double-layer potentials (which we shall sometimes refer to as SLP and DLP) satisfy 
the equation Cau = 0 with zero Cauchy data. However to solve the DNR problems, one needs 
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to study the boundary traces of these operators. To this end, one considers the boundary singular 
integral (or BSI for short) 


/C(/)(P,t) := lim/C,(/)(P,<), P&dn,t>0, 


where 


:= [ [ d^r{P,Q,t- s)f{Q,s)dQds 

Jo Jan 


10 Jan 

n pt-e 

= 51 / I a^j{Q)nj{Q)dyT{P,Y,t-s)i^^^f{Q,s)dQds 




10 Jan 


=:J2icy^{f)iP,t), £>0. 

Remark 4.1. Note that one uses the prineipal value in the integral defining the boundary singular 
integral beeause the points P and Q in the integrand are both on the boundary and can get very close 
to each other resulting in an undesired behaviour in the exponential function hidden in the integral 
kernel of 1C, when t and s are close to each other. The principal value is not needed in the integral 
formulas for the single and double-layer potentials because the point X is an interior point while 


Q is on the boundary, hence they are separated. Also, as we shall see in Proposition 4.6 below, it 
makes no difference if we consider the principal value in “time” or in “space”. 


Remark 4.2. We also need to consider the adjoint operator 


pt — € P 

IC*{f){P,t) := / / aij{P)nffP)d,,T{X,Q,t-s)\^^^f{Q,s)dQds. 

.-I Jo Jan 


i,j=i ' 


Note this presentation is valid thanks to our assumption A* = A. All the results that we are going to 
prove for K. are also valid for 1C*, due to the same behaviour of their corresponding integral kernels. 

4.1. boundedness of BSI. For the application of Fredholm theory in showing the invertibility 
of the relevant boundary integral operators, the following boundedness result is crucial. 


Theorem 4.3. Assume that 


(41) 


(44) 


hold. Let e > 0. 


Then, 


ll^e(/)||L^(S^) < ||/||l=(s^), fCL \ S ^). 

Proof. The idea behind the proof is as follows. First, one takes the Fourier transform in time of 
fCeif) and rewrites the resulting operator as an elliptic boundary singular integral plus some error 
terms. For the elliptic part which contains cancellations, we take advantage of the results in the 
elliptic theory, previously established in |S], while the error terms will be controlled by the Hardy- 
Littlewood maximal function Aidn- Finally an application of Plancherel’s identity allows us to 
return to the original operator. Now, let / G Lf{Soo) and i, j = 1,... ,n. For every P G dCl we have 


^oo nt — e 


(18) /CF(/)(P,r) = 


>0 Jo Jan 


a^ffQ)nj{Q)^y.T{P,Y,t- s)|y,^Q/(Q, s)e dQdsdt 


Ian 


fiQ)nj{Q) / f{Q,s) / dyT{P,Y,t- s)\^^^e '■'^*dt dsdQ 


' s+e 


= / a^ffQ)nj{Q)f{Q,T) / dyT{P,Y,C)\^^^e 

J on ^ J e 


dQ. 


Then we split the above integral as follows: 


■■= hif) + W) + W) + hif) + hif), 
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where 


lann{V^<\p-Q\<^} 


ann{|p-Q|<v^} 


a,j{Q)nj{Q)f{Q,T)dyT{P,Y)\^^^ dQ, 

pOO 

aij{Q)nj{Q)f{Q,T) j dyV{P,Y,C)\^^Qe~"^dC, dQ, 

dQ, 


I a^j{Q)nj{Q)f{Q,T) / dyT{P,Y,C)i^^^e dC, 

-'ann{|p-Q|>v^} '-Jo 

f a,,{Q)n,iQ)fiQ,r)\ a,T(P, y, dC 

5nn{|P—Q|>max{v/£,—^}} ^Jo 


l2{f){P,r) = 
h{f)iP,r) = - 
h{f){P,r) = 

h{f){P,r) = 

First of all, observe that 

Mf)iP,T)\<2V’^{f{-,T)){P), 

where /C*’-' represents the elliptic boundary singular integral given by 


dQ, 


9nn{V^<|p-Q|<^} 


a^J(Q)nj(Q)f{Q,T) I 9j^T(P, F, (e - 1) dC dQ. 




T<P’^{9){P) ■= sup 
<5>0 


/ aij{Q)nj{Q)dyY{P,Y)\^^Qg{Q) dQ . 

/onn{|p-Q|>5} 

Thus, the L^(S'oo)-boundedness of the integral di follows from [HI Theorem 3.1]. Next, we deal with 
the remaining integrals I 2 to 
From Lemma [3.31 it follows that 

dC 1 


^ 9,T(P,y,C)|,^«e—«dc|<^ 


Hence |(T3)| yields 


l^2(/)(P,r)| 


< 


e(" i)/2 dann{|p_Q|<^} 


Q{n+l)/2 s(n-l)/2- 

\KQ,T)\dQ<MoM;T)){P). 


Lemma |3.3| once again yields 


a,T(p,y,c)|,^„e—^dc 


< 


< 


re g-c|P-Q|Vc 

L 

p-c\P-Q\Qe 


^dC 


1 


Jc\P-Q\ye 




\P-Q\ 


n—1 


Next, we apply Lemma 2.1 to obtain 

|/3(/)(P,r)|< 


/ann{|p-Q|>v^} ^ 


-c|P-Q|Ve _ 

\fiQ,r)\dQ<Manifi;r)m. 


On the other hand, using Lemma 3.8 (in) for any N G N \ {0}, and using Lemma 2.1 it follows 
that 


1 


|/4(/)(P,r)|< . 

rl dann{|p-Q|>^=} \P - Q\ 

Vkl 


dQ < Manifi; r))(P). 


Finally, we once again use the fact that |e — 1| < I^CI^j for all 0 < /I < 1, and therefore Lemma 
|3.3| yields 

(19) 


r dyY{P,Y,(:)i,^^{e-^^^-l)dC <\Tf [ 

Jo Jo 


°° e-c|P-Q|VC 


^(n+l-2/3)/2 


dC 


-ld 


\p-Q\n-l- 2 p 

Using this estimate and Lemma [2.1] we obtain 


g-s^(n-3-2/3)/2^^ 


-1^ 


|P-Q|"-i-2d- 


\h{f){P,r)\<r^ [ 
Ja 


< Manifi;r)m. 

ann{\p-Q\<^} \P - Q\ ^ 
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Summing all the pieces together, the boundedness of the Hardy-Littlewood maximal function 
and Plancherel’s theorem yield the desired result. □ 

Remark 4.4. The Banach-Steinhaus theorem, together with the above uniform Lf boundedness 
and the pointwise convergence of IC^f for functions in C'“(S'oo), yield the convergence of in 
the norm to the bounded operator 1C. Furthermore, based on this fact and on the estimates 
forT{P,Q,t), which are of the same nature as in the constant coefficient case, standard Calderon- 
Zygmund theory yields that the operator JC is bounded on L^(S'oo) for any 1 < p < oo. 

As a consequence, we obtain the following pointwise convergence result: 

Corollary 4.5. Assume that \{Al)\ - (A4)| hold. The operator given by 

/C(/)(P,t) :=sup|/C,(/)(P,t)|, Ped^},t>o, 

s>0 

is bounded in L^(5'oo). Hence, for every f G L^(5'oo) the limit 

^(/)(-P>0 := lim P G dn, t > 0, 

e-fO 

exists almost everywhere. 

Proof. Define lf{f){P,t) as 

J j Q ^X{t—s>e} ('S) X||p_Q|2^t_s>£| (Qj dQds, 

and let I^{f){P,t) := KPd{f){P,t) — I^{f){P,t). For the sake of simplicity, here and in the rest of 
the proof, we suppress the dependency on the i,j parameters. Then, using the estimates for clyF 
from Lemma |3.3| and following a similar argument as in the proof of Theorem |4.3| yield 


( 20 ) 


sup \llif)iP,t)\ < Man {M,if){t)) (P). 


e:>0 


Now observe that for each {Pi,ti) G S{P,t,e) where 

S{P, t, e) := {b{P, \/s/S) n x (t — e/3, t + e/3), 

we have 


( 21 ) 

=^(/)(-Pi,ti) -/C (/(Q,s)x||p_q|2_^j_^<^|(Q,s)) (Pi,ti)+ 

+ j s) - K{Pi,Q,ti - s))x.{ip_Qi2^t_,.^,^j{Q,s)f{Q,s)dQds, 

where K{P,Q,t) = aij{Q)nj{Q)dyT{P,Q,t). 

We claim that, for any 1 < q < 2 

( 22 ) 

sup |4'(/)(P,t)| < Man{Mi{JCf){t)){P) + (Mon(A^i(|/r)(t))(P))'/'' + A^i(Xan(/)(t))(P). 

e>0 

Averaging the second term on the right hand side of pT| over S{P,t,e), Holder’s inequality with 
1 < q < 2, and the M boundedness of the operator K, (see Remark 4.4 above) yield 


|5(P,t,e)| 


-1 


'S{P,t,e) 


< 


/S{P,t,e) 


K. (/(Q,s)x||p_Q|2+t_,<e}((3,s)) (Pi,ti) 
^ (/(Qj s)x|lP_Q|2_|_t_s<£|(Q, s)^ (Pi,ti) 


dPi dti 


1/9 


dPi dti 


< 


\ 1/9 


/ \f{Q,s)\^ dQds] <{Man{Mi{\fit)n)iP)f'^, 

\ Qp+t— s<e} / 

with constants independent on e. Also, clearly we have 


\SiP,t,e)\-^ I \ICif){Pi,h)\ dPMh < ManiMiilCf){t)){P). 

JSiP,t,s) 
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Thus it remains to show that the last term on the right hand side of (211 is bounded uniformly in 
e by Ali(Alan(/)(P))(<). To this end we observe that the aforementioned term is bounded by a 
constant times 

I I \dyT{P,Q,t- s) - dyT{Pi,Q,t- s)\\f{Q,s)\ dQds 


idn Jo 


pt-e+\P-Q\^ 


+ 


dy T{Pi,Q,t-s) - dy T{Pi,Q,ti - s)| |/(Q,s)| dQds. 


/an Jo 


Therefore, using the mean value theorem and Lemma |3.4| we have that the expression above is 
bounded by 

I- i-t-e+lP-Qf g_c|Z-Q|2/(t-s) 

sup !/(<?> 


/an Jo 




f.t-e+\P-Q\^ 


+ £ 


/an Jo 


g-clPi-QlV(r-s) 


where [a, b] denotes the segment connecting the points a and b. 

The analysis of the two terms above are quite similar, so we confine ourselves to deal with the first 
one. To this end, we decompose the boundary integral in the first term above into two integrals - 
one over cirin ||P — Qj < •\/£/2| and another over its complement. The first resulting integral can 
be bounded by 


JaQn{\p-Q\<^}Jo (t-s)("+2)/2 


- s )("+ 2)/2 


ds 


where the last step follows from Lemma 2.1 Using the triangle inequality we have that \Z — Q\ > 
C \P — Q\foT Z G [P, Pi]. Hence the term which is defined by integrating over |P — Q\ > y/e/2 
is bounded by 


/■ p-c\P-Q\‘/(t-s) 

Ve / 77-7771/ , 77 - in/ 2 - 1/4 

JO ~ aann{|p-Q|>y772} U ~ s) / / 


-clP-QlV(t-s) 


<5/i 


r r \f{Q,s)\ 

lo (t - yann{|p-Q|>v77^} |P - 

<Mi{Mon{f){P)m, 


dQ ds 


where the last step follows once again from Lemma |2.1| This ends the proof of the claim for 

iKf)iP,t)- 

Now , using (20) and (22), the boundedness of the Hardy-Littlewood maximal functions and Remark 
4.4 we obtain the desired boundedness result for /C, from which the pointwise converge follows at 


once. 


□ 


In connection to the jump relation for the double-layer potential, the following proposition will 
prove useful. 


Proposition 4.6. Let f G L'^{St)- Then for a. e. {P,t) G St, 


lim [ [ d,jr{P,Q,t — s)f{Q, s)dQds = lim [ [ 

JdQ Jd 

Proof. Since 

[ f d,,r{P,Q,t-s)f{Q,s)dQds- [ f 
Jo JdQ JO Jd 


0 f 2 n{|P—Q|>v^} 


9i.r(P, Q,t- s)f{Q, s) dQ ds. 


lo ./9nn{|p-Q|<v^} 

rt r 


0 ./aan{|p-Q|>v/^} 
d,,r{P,Q,t - s)f{Q,s)dQds 


d,yT{P, Q,t- s)f{Q, s) dQ ds 


It—e JdQ,r]{\P—Q\>^} 


d„T{P,Q,t- s)f{Q,s) dQ ds =: Ief{P,t) + JJ{P,t). 
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It will be enough to show that, for all / G C^{St) 


(23) 


|sup|4/|||l2(s„„) WfWL ^ s ^), 

s>0 


(24) 

and 

(25) 

Now 

(26) 


SUp|JJ|||i2(s_^) < ||/||l2(5^), 
£>0 


lim Isf = lim J^f = 0. 

S —^0 S —^0 


would follow, if we could show that 


nt — e 


sup| / / a^j{Q)nj{Q)dy.r{P,Y,t-s)\^^^f{Q,s)dQds\<Mi{Mo!:i{f){P)){t). 

e>o Jo Jann{\p-Q\<^} 


But then Lemma 3.3 yields that the left hand side of (26) is bounded by 

pt—e r 


lo JonnilP-Ql< 2 V^} {t - s)("+i)/2 ^ ^ _ ,)(„+i )/2 

<Mi{Mon{f){P)m, 


ds 


where in the last step we applied Lemma 2.1 This shows (23). To prove (|24|, it is enough to show 
that 


(27) 

sup| [ [ a^J{Q)nj{Q)^yT{P,Y,t-f{Q,s)dQds\<Mi{Manif){P))it). 

e>0 Jt — e daQn{\P—Q\>\/€} 


Hence lemmas 3.3 and 2.1 yield that the left hand side of (27) is dominated by 

rt p c\P-Q\^/{t-s) rt r 


lt-eJann{\p-Q\>y^} (t-s)("+b/2 


-\f{Q,s)\dQ ds 


< 


f f \fiQ,s)\ 

It-e Jann{\p-Q\>y^} 


dQ ds 


< - / Manifi; s))(P) ds < Xi(Xoa(/)(P))(<), 
£ Jt-e 


and (24) follows easily from this. 


Since the proofs of (25) are similar, we confine ourselves to the proof of lim£_>o 7e/(P, t) = 0. 
To this end, without loss of generality, we translate the limit from the point P to the origin and 
hence aim to prove that lime_>o/ e/( 0, t) = 0. Note that in dealing with this limit, we can locally 
write P = {P',(p{P')), Q = {Q' t^p{Q')), P',Q' G for a certain Lipschitz function ip with 

\^{P')\ < \P'\^{\P'\) where w > 0, ||a;||Loo < ||V(/3 ||loo < 1 and limj_,.o+a;(t) = 0. Then, we have 
that 


4/(0, t)= ^ 


i,i=l ■ 


in{|Q'P + |v(Q')P<e} 


9^T((0,0),y,t-s)| 


V = (Q'.<e(Q')) 


,,,P*j (Q',s) dQ' ds, 


where Fij{Q',s) := aij{Q',ip{Q'))nj{Q',(p{Q'))f{Q', (p{Q'), s) a/1 + \V(p{Q')\'^. Observe that, since 
Fij G LF'{St), in order to show that limE_,.o 4/(0, t) = 0, it is enough, by a standard density argu¬ 
ment, to show the result for Fij{Q' , s) = g{Q') h{s) where g and h are smooth compactly supported 
functions. Moreover, using Lemma 4.4] we have that limE_>o/e/(0, t) = lim£_>o 7e/(0, t) where 


(28) 4/(0,4:=^ 


^ pt-e 




ayT((0,0), y, t - s)|^^ , ,q,„5(Q') Ks) dQ' ds. 
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Now we split (28) into the following three pieces 

pt—s r 


fO 




rt—e 


(a,T((0,0), F, t - - 5,t(( 0, 0), r, t - 

X g{Q') h{s) dQ' ds 

%r((0, 0), Y, t - {giQ') h{s) - 5(0) h{s)^ dQ' ds 


Jo JR"-in{|Q'|<V^} 

+ / 9{0)Hs) [ dyV{{Q,{)),Y,t- s)\ dQ' ds 

Observe that the mean value theorem and Lemma |3.3| yield that 

rt—e 


JO JR-' 


e-c|QT/(i-s) 

77- \(n+2)/2 


< 


Is JR"-in{|Q'|<v^} 

<^/i 


IQV(IQ'I) 


e-dQT/s 

s ("+ 2)/2 


dQ' ds 


'-in{|Q'|<Vi} 


u(VilQ'l)dQ'){j 


-in{|Q'|<i} 

Using again Lemma |3.3| and the mean value theorem to g, we get 


giu+2)/2 


ds —>■ 0, e —>■ 0. 


i/:y(ou)i< 


pt—e 


-c\Q'\Q(t-s) 


rt „-c\Q'\Qs 


< 


y-in{|Q'|<Vi} (t-s)(”+b/2 

dQ' 

in{|Q'|<V?} '■ 


^n{\Q'\<V^}Js s' 


nl2 


ds dQ' 


|Q'|"-2 


e-V"/2-2 dr) - Ve ^ 0, e ^ 0. 


Finally, note that r(X, Y,t — s) = ro(X, Y,t — s) + ri(X, Y,t — s), where 

^_(A-i(Y)(A-Y).A-Y>/4(t-s) 


To{X,Y,t-s) = Cr,^ 


{t — s)"/^(det 


See e.g. jHI for the details. Now we claim that 


/ %ro((0,0),F,t-s)| dQ' = 0. 

jR"-in{|Q'|<v^} 

This follows by using the same reasoning as in the proof of Lemma |3.3| and the oddness of the 
resulting kernel. 

Now for r 1 one has the estimate 


(29) 


=-c|A-y|V(t-s) 


|avri(x, ut - s)| < X(0,^)(t - s), 


where 0 < a < 1 is the Holder exponent appearing in assumption (A4). The estimate (291 follows 


from those in [5J p. 377], and once again from the same reasoning as in the proof of Lemma 3.3 
Therefore 


pt — e 


lim 
6 —^0 


T{|Q'|<V^} 


|a,Ti((o,o),r,t-s)|_, jdQ'dt = o 


since (291 allows one to apply the Lebesgue dominated convergence theorem. This in turn yields 


that lime_,.o/*'^ 3 /( 0 , t) = 0, and summing up, we obtain lime_,.o/£/(0, t) = 0 which concludes the 
proof. 

□ 
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4.2. boundedness maximal DLP. The estimates for non-tangential maximal functions of the 
layer potentials are crucial for establishing almost everywhere convergence of the solutions to the 
initial data as well as the jump relations, which will be used in the analysis of the invertibility 
problems. In analogy with the usual heat equation, the following estimate holds. 


Theorem 4.7. Assume that {Al) - (A4) hold. Then for f € Lf{Sao) one has 


mf))i 


< 


L^(S^ 


II/IIl 2 (Soo)> 


where {-jj denotes the non-tangential maximalfunction defined in 

Proof. Fix i,j = 1,..., n. Let P G d^l, X G 'y±(P) and set e :=\X — Pp. Then we can write 

= ici’^if){p,t) + Jiif)ix,p,t) + .h{f){x,Pfi) + Mf)ix,Pfi) + J4{f)ix,p,t), 

where 


Mf){X,P,t):= f " f 
Jo Jd 


aij{Q)nfiQ)f{Q,s) 


0 J dnn{\p-Q\<2./^} 

dyT{x, y, t - - 5,T(P, y, t - 


dQ ds, 


J2if)iX,P,t):= f " f 
Jo Jd 


aijiQ)nj{Q)f{Q,s) 


0 J dnn{\p-Q\>2^} 

dyX{X,Y,t- - dyY[P,Y,t- s)|^^ J dQds, 


Mf)iX,P,t):= f [ 
J t—£ J d 

Mf){x,p,t)-.= [ [ 

Jt-e Jd 


t-e Jd0.n{\P-Q\<2^} 
t 


>t—e 5r2n{|P—Q|>2y^} 
By Lemma |3. 3 1 we easily obtain 


aij{Q)nj{Q)f{Q,s)dyY{X,Y,t - dQ ds, 

aij{Q)nj{Q)f{Q,s)dyT{X,Y,t- s)|y^Q dQ ds. 


\Mf){x,p,t)\< f " f 

Jo Jd 


\f{Q,s)\ 

0 Jann{\p-Q\<2^} {t — 

<M4iMoUf){P)m, 


dQds<e^-^y^ r ^oUfi;s)){P) 

^ ^ Jo (t-s)("+l)/2 


where in the last step we applied Lemma 2.1 


For J 2 we first use the mean value theorem and then Lemma 3.4 (i), to get 
(30) \J2{f){X,P,t)\ 


rt—e 


< 


/ / \fiQ,s)\\x-p\ sup 

/o Jdnn{\p-Q\>2^} zeix,p] 


dQ ds 


Jo Jdi 


lo Jdnn{\p-Q\> 2 ^} ze[x,p]\Z - Q\^ i/2(-^_g)5/4 

where [X, P] denotes the line segment which connects the points X and P. Moreover, since |X—P| = 
^/e, if |P — Ql > 2^/e then jX — P| < |P — <51/2 and therefore 

\P-Q\ 

2 


\fiQ,s)\ 


dQ ds, 


(31) 


|y-Q| > |p-Q|-|y-p| > 


for \P-Q\>2^/E, Z G[X,P]. 


Hence, the last integral in (301 can be controlled by 

ot — S 


^ r |/(Q, 5 )| 1/4 M9a{fi;s)){P) 

Jo it-sf/^Jenn{\p-Q\> 2 V-e}\P-Q\^-^/^ ^ Jo (f-s)5/4 


ds 


<MfiMan{f){P)){t), 


where we have used Lemma |2. 1 1 twice. 

For J 3 , estimate (0 and Lemmas 3.3 and 2.1 yield 

ft „-ce/(t-s) 


/ t — ce/(t — s ) 

^—^^-^Man{f{; s)){P) ds < MfiMan{fm){t). 
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Finally, (31), Lemmas 3.3 and 2.1 yield 

nt 


\Uf){x,p,t)\< ( f 

Jt-e Jd 


a-c\P-Q\‘^/{t-s) 




■\f{Q,s)\dQ ds 


< 


f I 

Jt-e Jdi 


\f{Q,s)\ 

|P-Qh+i 


dQ ds 


<1 


t—s J dQr\{\P—Q\>2y/e} 

* Manifi; s))T) ds < MiiMon{f){P)m. 


In conclusion, we obtain the pointwise estimate 

sup \V^’^f){X,t)\<r’^{f){P,t) + Mi{Man{f){P)m, {P,t) e S^. 
xej±{P) 


Hence to end the proof of this proposition, we use Corollary 4.5 and the L^-boundedness of the 
Hardy-Littlewood maximal operators. □ 


Remark 4.8. Since Lemma 3.3 part (i) and Lemma 3.4 part (i) yield the same estimate (as far 
as the decay in \X — Y\ is concerned) for the second derivative of d^T(X,Y,t) and i9^yr(X, L, t), 
using a similar argument as in the proof o/Theorem |4.7[ one has 


(32) 


(V5(/))^ 


< 






This will be important for the invertibility of the BSI associated to parabolic equations. 

4.3. The jump relations. The discontinuity of the double-layer potential in the non-tangential 
direction across the boundary is reflected in a precise way in the following jump relation. 


Proposition 4.9. Assume that (Al) - (44) hold. Let f S L^{Sao)- Then, 

lim V(f)(X,t) = ^]^f(P,t)YK.{f){P,t), a.e. {P,t) G Soo- 

X^P z 

xe'r±{P) 

Proof. Taking into account Proposition |4.7| and using standard techniques, it is enough to see that, 
for every / G L^(S'oo) such that / G C(P(Sac,), 

hm ^){X,T)=Tlmr)+]^){P,T), a.e. PG^fl. 

x^p z 

XG7±(P) 

Let i, j = 1, ... ,n, P G dLl and X G 7 ±(P). As in ([T^, we can write 

V^){X,t)= [ a,j{Q)n,(Q)f(Q,T)\ r dyYiX,Y,0^,^^e-^^^d(:]dQ 
Jan '-Jo 


Ian 


aij(Q)nj(Q)f{Q,T)dyT{X,Y)\^^Q dQ 


+ a,yiQ)n,(Q)f(Q,T) / a,T(A,y,C)|,.Q(e-*"« - l)dC dQ 

Jan '-Jo 

=:A^)\f)(X,r)+A^\f)iX,r). 

Recall that, from the elliptic jump relation (see for example [HI Theorem 4.6]), 


1 ^ /* 

(33) lim A1’^(/)(A,r) = T x/(C t) -f liin / 

X—>P ^' Z £—>-0 a 

XG7±(-P)»d=l ° 


ann{|p-Q|>v^} 


d,,r(P,Q)f(Q,T)dQ. 


Moreover, 

(34) 


hm Ay^(f)(X,T)=Ar{f)(P,T), 

X—>p 
XG7±(P) 


by the Lebesgue’s dominated convergence theorem. Indeed, estimates (191 and (|^ yield 


\JCfHf)iX,T)\ 


< 


\fiQx)\ 


Ian \X - Q|"-i-2/5 


dQ< 


dQ 


I an 


|P- Q|"-l-2/3 


< oo. 


where it was important that / is bounded. Observe that the constants involved in the above 
inequalities might depend on r but not on X. 
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It remains to see that the second term in the right hand side of (331 plus the sum of the right hand 
side factors of (341 equals JC{f). Since the integrals A 2 ^ (f) are absolutely convergent, we easily get 


lim 

S —^0 


dnn{\P-Q\>y/^} 


= lim 
€—^0 


d,T{P, Q)f{Q, T)dQ+J2 ^ 2 'ifKP, r) 


/onn{|p-Q|>^/^} 
On the other hand, by Proposition |4.6| we have that 

rt 


d^r(F,Q,T)f(Q,T) dQ. 


/C(/)(P,t) = lim / f 


0 Jann{|p-Q|>v^} 

Hence, proceeding as in (18 1 , we conclude that 
^{f){P,T) = lini 


dyT{P, Q,t- s)f{Q, s) dQ ds. 


d^T{P,Q,T)f{Q,T) dQ. 


Jann{|p-Q|>v^} 

Notice that we have used (see Corollary |4.5| ) 

/C,(/)(P,-) ^/C(/)(P,-), inp2(0,oo), 

and the continuity of the Fourier transform in P^(0, 00 ). 

For the normal derivative of the single layer potential one has the following jump relation. 


□ 


Proposition 4.10. Assume that (41) - (44) hold and let f G L^(S'oo)- Then we have 
a,(5(/))±(P,t)=±i/(P,t)+/C*(/)(P,t), u.e. (P,t)e5oo, 
where is defined as in (|^ and (•)^ is defined as in ( |10[ ) and ©• 

Proof. We observe that the integral kernel of d,yS is similar to that of DLP’s and hence verifies the 
same estimates. Therefore proceeding as in the proof of Proposition |4.9| we will obtain the desired 
jump relation. □ 

As a consequence of the jump relation above, and the definition of Vt one sees just as in the case 
of constant coefficient operators that 

Vt(5/)+=Vt(5/)-, 
where Vt is defined as in (|4]), see also [2]. 

Furthermore, taking the limit in the integral defining the SLP, using the estimate (3.1) and the 
Lebesgue dominated convergence theorem yield that 

{Sf)+ = (Sf)-. 

4.4. boundedness of the maximal fractional time-derivative of SLP. In connection to 
the problem of invertibility of the single-layer potential and the regularity problem, we would also 
need an estimate for the fractional derivative of the single-layer potential. The fractional time- 
derivative of the single -layer potential is defined as 

DP^Sif)iX,t):= f [ D]^^T{X,Q,t-s)f{Q,s)dQds, X G n, t > 0. 

Jo Jan 

Our goal is to prove the analogues of the estimates derived for the DLP, for the fractional derivative 
of SLP defined above. To this end we need a series of results showing the boundedness of various 
operators related to Let f 1 —>• K be Lipschitz and define a kernel JC by 

jr(p', Q', t) := py^r((p',<p(p')), (Q', ^m),t). 

Since ip is Lipschitz it follows that Jtf{P',Q',t) satisfies all the estimates in Lemma |3.5| if we 
replace \X — Y\ in the right hand sides of those estimates by \P' — Q'\. Now let r{P',Q') := 
■\/\P' — Q'p -I- |<^(P') — </5(Q')P for any e > 0 define for / G x (0, c»)) the operator 

Te, acting on / by 


TJiP',t)= f [ JY{P',Q\t 

Jo Jr{P',Q')>e 


s)f{Q',s)dQ' ds. 
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One now proceeds by investigating the Tg{f) and Tf := lime_>.o Ts{f) for all / G x (0, oo)), 

however in the proposition below, we view those / whose domain have been extended to x K 

in such a way that f{Q', t) = 0, for t < 0. This enables us to take the Fourier transform of f{Q', t) 

in t, which will prove useful. 

Proposition 4.11. Let K" ^ := K" ^ x (0,oo) and f G ^). Then one has 

(i) ^ II/IIl 2 (][{"- 1 ) where the constant hidden in the right hand side of this esti¬ 

mate doesn’t depend on e. 

(ii) The limit lime_>oT e(/) exits in the norm and therefore it is legitimate to define 

Tf as this limit. 


Proof. To prove (z), using (15), we observe that the time-Fourier transform of T^f is bounded by 


\Tef{P',r)\< 
Now Lemma [3 . 7| yields 


'r(P',Q')>e 


\^/^\f{{P',ip{P')),{Q',g,{Q')),T)\\f{Q',T)\dQ\ 


\f{{P',p{P')),{Q\<f{Q')),T)\< 


min{l,(|r|r(P',Q')')-'^} 


and therefore we can proceed as in the proof of Lemma |3.9| with q = 1/2 to show that 


/ |r|'/2|r((P', g,{P')), iQ\ g,{Q')l r)| dQ' < 1. 

Thus, Lemma [2.2| and Plancherel’s theorem imply the boundedness of T^f. 

The proof of (zz) is standard since it amounts to show that the sequence {?//} is Cauchy in 
This is done by using Fubini’s, Plancherel’s and Lebesgue’s dominated convergence 
theorems. The details are left to the reader. □ 


Now if one uses the fact that for {P',t) 9 supp (/) one has 


Tf{P',t) = 


JiT {P', Q' ,t — s) f{Q', s) dQ' ds, 


then thanks to Lemma 3.5 one can follow the exact same steps in |21 Appendix A, Proposition 
A.3] to show that |jr/||^ 2 (R--i) < ||/|Il 2 (r^- 1 ). Moreover, setting T*/ := smp^^oTs{F){P',t) and 
following the same argument as in [3 Appendix A, Theorem A.5], one can show that 


(35) 


\n 




Il2(r" 


We would like to emphasise once again that these estimates follow using the same method as 
in the constant coefficient case in |3, because the proofs there are all only dependent on estimates 
for the kernel JF{P',Q',t — s), which due to Lemma 3.5 are the same as those in the constant 
coefficient case. As usual the maximal function estimate (35) yields that, for / G one has 


Tf{P',t) = \im [ [ J(f{P',Q',t-s)f{Q',s)dQ'ds, 

Jr{P',Q')>s 

for almost all {P' ,t) G 

Now using the definition 0 of the boundary of the Lipschitz domain LI, one can transfer the 
operator T to Sao, by setting 


T/(P,f) = lim f [ JF{P,Q,t — s)f{Q',s)dQ'ds. 

Jdnn{\p-Q\>e} 

Therefore using the boundary coordinates {P', ip{P')), one sees that 

ffiiP', g,{P')),t) = T{f Vl + |V^|2)(P', t), 
and hence all the results above for T and T» are also valid for T and T*. 

Finally, following the same arguments in [21 Propositions A.7 and A.8], line by line, one has the 
following two results: 
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Proposition 4.12. Assume that\{Al) - (44)1 hold. Then 




< 


LHS^) 


WfhHS^), f&L^{Soo). 


To establish the jump relation for the fractional derivative of SLP, one also needs the following 
result: 


Proposition 4.13. Assume that (41) - (44) hold and let f G L^{Sao)- Then for a.e. {P,t) G Sac, 


lim Dy'^S{f){X,t) =Tf{P,t) = \im f f Dl^^r{P,Q,t - s)f{Q, s) dQ ds. 

Jdnn{\P-Q\>V^} 


x^p 

A:e7±(P) 


The proof of this proposition is once again the same as the corresponding result in the constant 
coefficient case in Appendix A of j^. From Proposition 443[ it follows at once that 


{Dy^Sf)+ = {Dy^SfY 


For the invertibility of the SLP discussed in Subsection |7.2| the following boundedness estimate 
plays a crucial role. 


Theorem 4.14. Assume that Ai and A 2 satisfy the conditions |(41) - (44)| and that Ai = A 2 on 
dTl. Then one has for 0 < T < 00 

(36) ||V(5a, - SaYIWlhs^) + - SaYIUhs^) + IK^a, - SaYIUhSt) 

<||4i-42||V2||/|U.(s,). 

Proof. Since 

Dl^\SA,-SAYfiX,t)= f [ Dl/^{TAAX,Q,t-s)-TAAX,Q,t-s))fiQ,s)dQ 

Jo Jan 


ds, 


using Lemma 3.6 we observe that the Schwartz kernel of the integral operator dI^^(Sai — SA 2 ) is 
given by 

^ nOO p 

KAX,Q,t-s) = Y] / / d,TAAX,Z,u)d,^DYTA 2 iZ,Q,t-s-u){AAZ)-A 2 iZ))dZdu. 

ij=l "'0 

Now we observe that since Ai = A 2 on did for any point Q G 11 and any Z G fl we have 


AAZ) - A2(Z) = AAZ) - 4i(Q) + 4i(P) - A2(Q) + 42 (g) - 42(Z) 

= Ai(Z) — 4i(g) + 42(Q) — A2(Z). 

Hence, using the fact that 4i and 42 are Holder-continuous matrices, we have that for any point 
Q G n and any Z G H, one has that 


(37) 


\AAZ)-A 2 iZ)\<\Z-QY, 0<a<l. 


Taking this fact into account. Lemma |3. 3 [ and Lemma 3.5 part (ii) yield that 
\KAX,Q,t-s)\ 


t-s r g-clA-ZiVij 


^ — c.\Z—Q\^l(t—s—u) 


<U^-A2fJY 

<||4i-42||L/^ [ 

Jo 

<||4i- 42||L/" r f 

Jo [t-s- uy/^ V7|z|>(t_s-„)i/2 


0 Jk^ u("+ 1)/2 (t_s_„)3/2|^_Q|n-l-a/2 

_ 1 

g u("+l)/2 [t — g — uy/"^ 


dZ du 


g-c|X-Q-Z|V« ^-c\Zj/{t-s-u) 


dZ I du 


|^|n-l-a/2 

^-c\X-Q-Zj/u ^-c\Zjjit-s-u) 
| 2 '|"-l-a /2 


dZ 1 du 


g-c|X-Q-Z|V« g_c|Z|V(t-s-u) 


+ \\A,-A2\\Yf 7.- dZ]du 


i+l)/2 (t-s- uy/^ \ J\z\<{t-s-uy/^ 


| 2 '|"-i -“/2 


:=I + J. 
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Now using Chapman-Kolmogorov formula (12) we have 

rt—s 


= \\A,-A,\\l/^(^J^' 


du 


^-c\X-Q'd' jt-s 


g y(ra+l)/2-ra/2^^ _ g _ y^3/2+(n-l-a/2)/2-n/2y _ g'jn/2 

t-s , ^-c\X-Q-Z\^/t-s 


l/2(i_5_y)l-a/4; (t-s)"/2 


= Bil/2,a/mAi-A,\\]^^ 


u 

du N ^-c\X-Q-Z\^jt-s 

yl/2(l _ u)l-«/4 J _ g)(«+l)/2-a/4 
g-c|X-Q|Vt-s 


(t _ s)(n+l)/2-a/4’ 


where i?(-, •) is Euler’s Beta function. 
To estimate J we observe that 


r ^-c\X-Q-Zf ju ^-c\Z\^Kt-s-u) 

/|Z|<(t-s-ii)i/2 |Z'|’^-l-“/2 


dZ 


-E / 

J=0 “'2-3-l(t-S-«)l/2<|2|<2-J(t-S-«)l/2 


g-c|A-Q-Z|V« g-c|Z|Vh-s-ii) 


|^|n-l-a/2 


dZ 


^ ^ 2 ^(n- 2 ) g-c2-V2(i _ 5 _ ^)-(n-l-a/2)/2 

j=0 


-c\X-Q-Zf/2u ^-c\Z\^/2(t-s-u) 


/2-2-l(t-S-tt)l/2<|Z|<2-3(t-S-«)l/2 


< _ g _ y)-("-l-a/2)/2 f ^-c\X-Q-Zf/2u ^-c\Z\^/2(t-s-u) 

JB" 


Therefore, using Chapman-Kolmogorov formula (12 1 again we obtain 

du 


J<I|24i-A2||^^(^ 


2+1)/2-"/2(^ - S - u)3/2+(n-l-a/2)/2-n/2y (2t - 2s)"/2 


,-c|X-Q|V2(i-s) 


g-c|Y-Q-Z|V2(i-s) 
( 2 t- 2 s )"/2 


< 


i?(l/2,a/4)||Ai-Al2||^/^ 


„-c|X-Q|V 2 (t-s) 


Thus 


\KiiX,Q,t-s)\< 


(i_s)(n+l)/2-a/4- 

Pi -242 II 


(|X-Q|2-Ht-s)("+l)/2-a/4- 

For the operator V( 52 ii — Pa)! observe that its Schwartz kernel is given by 

^ pOO p 

K 2 {X,Q,t-s)= V / / d,,d,TA,{X,Z,u)d,^rAAZ.QA-s-u){Ai{Z)-A 2 {Z))dZdu. 


This, Lemma 3.4 part (i). Lemma 3.3 and (371 yield that 

|K2(pg,t-s)| < [ [ \d,,d,TAAX,Z,u)d,^rAAZ,Q,t-s-u)\\A,{Z)-A2{Z)\dZdu 

Jo VB" 

■ dZ du 


< Pi-242^2 r f 

Jo Jw 


{\Z - X|2 -H u)("+ 2)/2 (\Z-Q\2+t-S- u)(’^+l)/ 2 -a /2 


< 


Pl-PIlL^' 


(p-Q| 2 -Ht-s)("+l-“)/ 2 ' 

Finally for — p^ we observe that its Schwartz kernel is 


^ pOO p 

K3(X,Q,t-s)= V / / d,TAAX,Z,u)d,XA,{Z,Q,t-s-u){Ai{Z)-A2{Z))dZdu. 
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This, Lemma 3.3 and (371 yield 


\K^{X,Q,t-s)\< [ [ \d,TAAX,Z,u)d,XAAZ,Q,t-s-u)\\Ai{Z)-A 2 iZ)\dZdu 

Jo Js.^ 

r f 

Jo Jr 


< 


,0 JR^ {\z - X|2 + u)(-+l)/2 (|Z - Q|2 + t - s - u)("+l)/2-«/2 

\\a,-A2\\]1^ \\a^-A2\\]L^ 


dZ du 


{\X -Q\^+t- s)("-“)/2 {\X - Q|2 + i - s)(n+l)/ 2 -(l+a )/2 ' 

In conclusion, combining the above estimates we get 


\Ki{X, Q,t-s)\ + \K 2 {X, QA-s)\ + \K 3 {X, Q, t 


^ \\M - A 2 \\]i^X{o,t){s) 

"" (|X - Q|2 + t - s)(n+l)/2-a/4 ’ 


Now observe that 


X(o.t)(s) ^ X(o.t)(s) 1 

(|X - Q|2 + t _ 5 )(n+l)/ 2 -a /4 - _ 5 )l-a /8 |p _ Q| („-l)/2-a/8 ' 

Therefore, since the estimate in the proof of Lemma |3.9| shows that 


sup 

P^dVL 


sup 

Qeao 


^ Ian 

dQ N 

|p_ Q|(n-l)/ 2 -a/ 8 j 

^ Ian 

dP \ 

|p_g|(n-l)/ 2 -a/ 8 j 


< 1 
< 1 , 


and since 


sup 

0<i<T 


sup 

0<S<T 



X(0.t)(s) 

(t_s)l-a /8 

X(0,t)(s) 

(t-s)l -“/8 


ds^ 

dtj 


sup 

0<t<T 


sup 

0<S<T 


ds 




< sup < 1 

0<t<T 


dt 


(i_s)l-a /8 


< sup (r 

0<s<T 


< 1 , 


Schur’s Lemma implies (36). 


□ 


5. Layer potential operators associated to the Fourier-transformed equation 

As mentioned earlier, the independence of the diffusion coefficient matrix of the time variable is an 
advantage here which enable us to perform a Fourier transform in the time variable in the parabolic 
equation and bring it to a parameter-dependent elliptic equation. However, one then needs to 
establish all the estimates we derived in the previous sections for the corresponding layer potentials 
of the Fourier-transformed equation. 

To this end, we define the Fourier-transformed single and double-layer potential operators by 

S"( 5 )(A) := [ f(A, Q, T)giQ) dQ, X e ft, 

Jan 


and 


D"(5)(A) 



d^f{X,Q,T)g{Q)dQ 



aij{Q)nj{Q)dyf{X,Y,T)\ 


y=q9{.Q) jq 


=:J2{wy’^ig)ix), Xen. 

i,j=^ 


We also consider the boundary singular integral 

















SOLVABILITY OF DNR PROBLEMS 


25 


where 


Ki9){P) ■■= 


/onn{|p-Q|>e} 


d^r{P,Q,T)g{Q) dQ 


= E 


a^j{Q)nj{Q)dyT{P,Y,T)i^^^g{Q)dQ 


ij^j^Jdnn{\p-Q\>e} 
n 

=: £> 0 . 

*j=i 

Remark 5.1. Note that we would also need to consider the adjoint operator, which is given by 


ri9)ip) -li-E 


ij^-^JdQn{\P-Q\>e} 


a^]{P)nj{P)^,,T{X,Q,T)l^^^g{Q)dQ. 


Once again, and (K”^) will satisfy the same Lf estimates. 


Remark 5.2. In sM&sec^zows |5.1| and |5.2| i/te constants hidden in the right hand sides of the estimates 
are all independent of parameter r. 

5.1. boundedness of the truncated Fourier-transformed BSI. The following theorem is 
one of the main tools in proving the invertibility of Fourier-transformed boundary singular integral. 

Proposition 5.3. Assume that \{Al)\ - (T4) hold and let e > 0. Then one has 

( 5 ) 11^2(90) < ||g||L 2 (aQ), g e L^{dVL). 

Proof. Let g G Lf{dTl) and i,j = 1,..., n. We write, for each P G dfl, 


mr’^(9)(p) = 


/ann{s<iP-Qi<^} 


aij(Q)nj(Q)g(Q)dyT(P,r)i^^^ dQ 


/ a^j{Q)ny{Q)g{Q) / dyT{P,Y,t)\^^^e dt 

f ^nn{|P—Q| >max{£, i— }} 


dQ 


/ann{6<|p-Q|< J=} 

\/1 t| 


aij{Q)nj{Q)g{Q) j dyT{P,Y,t)\^^^{e - 1) dt dQ 




=: h{g){P) + h{g){P) + l 3 {g){P). 


It is clear that 


\hig){P)\<2r’^{g){P). 


Moreover, following the same method as in the proof of Theorem 4.3 for I 4 and I 5 , we obtain the 
estimate 


\l2{9){P)\ + M9){P)\<M9n{9)iP)- 

Corollary 5.4. Assume that (Al)| - [(344) hold. The operator given by 

K"(g)(P) :=sup|K^( 5 )(P)|, P G dO, 

e>0 

is bounded in L^(5f2). Hence, for every g G the limit 

K^g)iP) := \im Kf{g){P), P G dll, 

£—>0 

exists almost everywhere, and it also defines a bounded operator in Lf(dTl). 

Proof. The proof is standard and goes along the same lines as that of Corollary |4.5 


□ 


□ 
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5.2. boundedness of the maximal Fourier-transformed DLP. For the Fourier-transformed 
double-layer potential one also needs to demonstrate the boundedness of non-tangential maximal 
function, in order to show the corresponding almost everywhere convergence and jump relations. 
To this end we have 


Proposition 5.5. Assume that \{Al)\ - (^4) hold. Then, 

(W{g))t gGL^dn). 

L‘^{dQ) 

Proof. Fix i, j = 1,... ,n. Let P G il, X G 'j± (P) and take e := |X — P|. We consider the following 
decomposition 

{Wy’^{g){X) = (KlJ’={g)iP) + Mg)iX,P) + Mg)iX,P), 

where 


Mg)iX,P):= / a,,(QK(Q)g(Q) 

./ann{|p-Q|>2£} 


dQ, 


J 2 {g){X,P) := / a^j{Q)nJ{Q)g{Q)^y.T{X,Y,T)\^^QdQ. 

J dGir\{\P-Q\<2e} 


For Ji we first apply the mean value theorem, secondly Lemma [3.8| (ii); next the relation (311 and 
finally Lemma [2.1| to get 


\Ji{g){X,P)\< \giQ)\\X-P\ sup 

J dnn{\p-Q\> 2 e} ze[x,p] 


< e / l5(Q)l sup ^ 

4ann{|p-Q|>2e} ze[x,p] \z - w\ 

<Moa[9){P). 


dQ < e 


/ann{|p-Q|>2e} 


|g(Q)l 

\P-Q\^ 


dQ 


To estimate J 2 we just need Lemma [3.8| (i), and estimate ([^, 

\J 2 {g){X,P)\< [ | 5 (Q)||a,T(X,y,r)|,^JdQ 

./ann{|p-Q|<2£} 

<e [ TY^^dQ<Manig)iP). 

JdQn{\P-Q\<2e} W - W\ 

Thus, 

sup |(D"y'^(g)(X)|<(Kf’"r’^-(5)(P)+A4aa(5)(^^), P & dn, 
AG7±(P) 

and the proof of this proposition follows from Corollary |5.4| 


□ 

Remark 5.6. Since Lemma |3.8| part (i) and {ii) yield the same estimate for the second deriv ative 
of dxT{X,Y,T) and i9^yr(X, Y, r), using a similar argument as in the proof 0 /Proposition 
one readily sees that 


5.5 


(38) 


(VS"(5))^ 


< 


L2(9n) 


llffl|L2(an), g G L‘^{dil). 


This will be important in the proof of the Rellich estimates in Section]^ 


Theorem 5.7. Assume that Ai and A 2 satisfy the conditions (yll)| - [{AA) and that Ai = A 2 on 
dLl. Then one has {uniformly in r) 


(39) 

and 

(40) 

for g G Lf{dLl). 


{{K\y-{K\J)g 

{{K\y* - {K\J)g 


< 


L^idd) 


11^1 - ^2||^^||g||L2(an), 


< 


L2(an) 


11^1 - ^2\yjy\\g\\L^(an), 


Proof. We only prove (401, since it will be the one we will be using later. Moreover (391 follows at 


once from ( |4(j| ). The kernel of (IK(^J* — (IK(^^)* is given by 

n 

^ {ay{P)ny{P)d,YAAX,Q,T)\,^, - al{P)ny{P)d,YAA{X,Q,T)\,^,). 
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Using the fact that a\j{P) = a^j{P), Lemma 
the aforementioned kernel is bounded by 


3.10 


Lemma 


3.8 


(ii) and (i), and (37), we see that 


n p 

\K{P,Q,T)\<\\A,-A2\\]i^ Y, / \d,,d,TAAX.Z,T)\,_,d,TAAZ,Q,T)\\A2{Z)-A^{Z)\dZ 






dZ 


< 


P 1 - 7 I 2 IIL/" 


\P - Z|" IQ - ^ |P - Q|"-i-“/2 ■ 

Therefore since supg^^Q fgn '^)l supp^^Q /an I^(^j Qj ''’)I both bounded by 

Schur’s Lemma implies the desired boundedness. □ 

5.3. The Fourier-transformed jump relations. Here we shall consider the jump relations that 
are valid for the Fourier-transformed layer potential operators. We start first with the jump relation 
for the Fourier-transformed DLP. 


Proposition 5.8. Assume that (Al) - {A4) hold. Let g € L^{dfl). Then, 

hm D"(5)(X) = t^5(^^)+K"(5)(^^), a.e. P € dn. 
x^p l 

X(^l±{P) 


Proof. From Proposition |5.5[ by standard techniques and density arguments it suffices to consider 
the case of 5 e Cf°{di}). 

Let i,j = l,...,n, P G dil, X £ j±(P) and r > 0. We have that 


y’^{g){X)= / a,,{Q)n,{Q)g{Q)dyT{X,Yy^^dQ 


IdQ 


+ / a^j{Q)nj(Q)g{Q) / 9yT(W U C)k=Q (e - l)rfC 


/an 


'-ao 


dQ 


=:A^\f)ig){X) + Ay{f){g){X). 
As in the proof of Proposition |4.9| we get 


" I 7 ~ 

hm ^ A^{\g){X) = T vff(P) + lim / d.r(P, Q)giQ) dQ. 

2 e ^0 7ann{|p-Q|>5} 


x^p 
X&1±(P) Li=l 


and 


hm A^\g){X)=A^\g){P). 
X&^±{P) 


Moreover, it is straightforward to see that 
lim / 


^^0 Jann{|p-Q|>£} 
which allows us to conclude this proof. 


5,F(P,Q)5(Q)dQ+ ^ A^i\g){P)=W{g){P), 

i j' = l 


□ 


For the normal derivative of the Fourier-transformed single-layer potential, we have the following 
jump relation which will be used in proof of the invertibility of boundary singular integrals. 


Proposition 5.9. Assume that (Al) - (A4) hold. Let g € Lf{dVt). Then 

a,(S"(g))±(P) = ±]^g{P) + {Wr(g){P), a.e. P G dn. 

Proof. Since the integral kernel of 9j/(S'^) is similar to that of the Fourier-transformed DLP and 
therefore verifies the same estimates, proceeding as in the proof of Proposition |5.8| we will obtain 
the desired jump relation. □ 

As a consequence of the jump relation above, and the definition of Vt one sees as before that 

VT(S"g)+ = Vt(S" 5 )-. 

Furthermore, taking the limit of the integral defining the Fourier-transformed SLP, using estimate 
(3.7) and the Lebesgue dominated convergence theorem we obtain 

{sy)+ = (sy)-. 
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6. Parabolic Rellich estimates 

6.1. Rellich estimates for the elliptic Fourier-transformed equation. 

Rellich estimates are the most basic tools in proving the invertibility of boundary singular integral 
operators appearing in the study of elliptic and parabolic DNR problems. In this section we pro¬ 
duce a systematic way to derive a family of Rellich-type estimates which will be enough for proving 
the invertibility of our boundary singular integrals. This is done by a rather general transference 
method which transfers estimates valid for operators with smooth diffusion coefficients to those for 
operators with Holder-continuous diffusion coefficients. 

We begin this section by recalling that Tq < oo will denote the diameter of and Qq is a cube that 
contains D. Moreover 2Qa will denote the cube with the same centre as Qa but with the double size. 

In what follows, it is important to note that, by multiplying the Fourier-transformed equation 0 
by u{X,t) and integrating by parts one obtains the following useful identity 


(41) 


Ida 


L{P,T)dy^w^{P,T)dP = ± [ i'^aij{X)diu{X,T)dju{X,T)+iT\u{X,T)\‘^\ dX. 

V / 


Now we start by proving some basic estimates concerning the solutions of the Fourier-transformed 
equations. It is important to note that in all the estimates that are established here (until the end 
of this subsection), the constants hidden in the right hand sides of the estimates are all independent 
of parameter r. 


Lemma 6.1. Let A he a matrix satisfying properties {Al) - (T4) Assume that Cau = 0 in LI. 
Then for every £ > 0 


(|Vu(y,r)p + |r||u(y,T)p) dY < (^e\d,,u^{P, t)\^ + 


2 , 


dP. 


Proof. The estimate follows easily by using (411 and Cauchy’s inequality ab < eaf /e, e > 0. □ 


Lemma 6.2. Let A be a matrix satisfying properties (Tl) - (T4) and assume that Cau = 0 in LI. 
Then 

[ \u{P,T)\^dP< [ (\Xu{X,T)\^ + \uiX,T)\^)dX. 

Jan Jn± ^ ^ 

Proof. Given a vector field H Integration by parts (i.e. divergence theorem) yields 

[ {H,N{P))\u{P,T)\^dP = ± [ \u{X,t)\^ {X-H)dX±Re'^ [ H^{X) dku{X, r) u{X, r) dX 

Jan Jn± ^ Jn± 

From this one deduces that 

(42) [ \u{P,T)\^dP< [ (|u(WT)||Vu(WT)| + |u(X,r)p)dX. 

Jan Jn± ^ ^ 


Therefore (42l and an application of Cauchy’s inequality to the term |u(X,r)||VM(X,T)| yield the 
desired estimate. □ 


Lemma 6.3. Let A G C°°(2Qn \ dLl) be a matrix verifying (13 1, then for each 6 > 0 one has 


[ |VH(X)||Vrt(X,T)pdX < [ 

Jn± Ja 


sup \Xu{X,t)\^ 
dfl (P) 


dP- 




\Vu{X,t)\‘^ dX, 


where 'ysiP) ■= 1+{P) C S LI : dist(X,911) < 5}, is a truncated cone with the vertex P G dLl, 
and ao is the parameter appearing in (131. 


Proof. We confine ourselves to the proof for Lt'^, as the one for LI is similar. Set 12+ := LI 
and write it as a union of the sets Ei and E 2 where E^ := {X G LI : dist(X, dLl) < 5} and 
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E 2 := {X G ft : dist(X, 9f2) > <5}. Then using (131 we have 


[ \VA{X)\\Vu{X,t)\^ dX = [ \VA{X)\\Vu{X,t)\^ dX + [ \VA{X)\\VuiX,T)\^dX 
*/ O J El J E 2 

^ f |Vw(X,r )|2 


dist(X, ^ j(i-Qo) 


/ E 2 


< 


\Xu{X,T)\^dX 
|Vu(X,t)|2 


< 


an 7{|Y-P|<(l+a)dist(X,0n); dist(A,an)<5} dist(X, “0 

f sup |VM(X,r)p( f dx \ ,^ 1 

aanxGTdP) 


dXdP- 


1 


,5(i-«o) 


|Vu(X,r)pdX 


/an xe-ysiP) 

<(5"-i+“« [ sup |Vu(X,r)pfiP + 
aan xe'jsiP) 


^X-Pl<s \X-P\^-o 
1 


dP 


5(l-ao) 


|Vu(X,r)|2dX 


□ 


For operators with smooth diffusion coefficients one has the following Rellich-type estimate. 


Proposition 6.4. Let A be as in Lemma 2.6 and assume that C^u = 0 in 12. Then, 


(43) [ (|Vn±(P,r)|2 + |r|KP,r)p)dP< / {\d,,u^{P,r)\^ + \u{P,t)\^) dP 

J 00 «/ 


and 

(44) 


[ |Vn±(P,T)|2dP< / (|Vrn(P,r)p + |r||n(P,r)|2 + 

Jon Jdn ^ 

+ ^n-i+«o f sup |Vu(X,T)pdP, 
Jon xgtPp) 




dP 


for any i5 > 0 and e G (0,eo(^))- 


Proof. First of all we observe that (41) yields 


(45) f (\Vu{X,t)\^ + \t\\u{X,t)\A dX < [ |5,,u±(P,r)||n(P,r)|dP. 

jn± ^ ^ Jdn 


Applying Cauchy’s inequality to the left hand side of (45) we also get 

(46) [ \T\-^\Xu{X,T)\\u{X,T)\dX< [ \d,,u^{P,T)\\u{P,r)\dP 

an± Jdn 


Now, given a vector field H and using the Fourier-transformed equation (17), an application of 
the divergence theorem yields 


n « 

(47) 5] / (iV(P),P(P))a,,(P)a,n±(P,r)a,u±(P,r)dP 

, JdVt 

n p 

'Y' / Hk{X) d^(dij{X)dku{X,T)dju{X,T)) -iTu{X,T)dku{X,T) 

AI — 1 -inL 


i j' = l ' 


= 2Re 




dX 


+ Y {X ■ H)dij{X)diu{X,T)dju{X,T) 
/n± 


i ,1 = 1 ' 


dX, 


where N denotes, as before, the unit normal vector to the boundary. Applying first the divergence 
theorem and then (45) to the first term on the right hand side of (47), we have 
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(48) 


n « 

|2Re ^ lHk{X)di(^dij{X)dku{X,T)dju{X,T)jdX\ 

n p 

= |2Re V / dMX)d,j{X)dku{X,T)dju{X,T)dX 

i j k—1 
n p 

+ Y / Hk{P)n^{P)dij{P)dju{P,T)dku{P,T)dP\ 

< f (^\d^^u^{P,T)\\Xu{P,T)\ + \d^^u^{P,T)\\u{P,T)\j dP. 


For the second and the third terms of (471, we use (461 and (451 and Cauchy’s inequality which 
yield 


(49) ir [ u{X,t) dku{X,T) dX + [ {X ■ H) dij{X) diu{X,T) dju{X,T) dX 

< (|a.,u±(P,r)|2 + |u(P,r)p + |r|i/2 |u| |a.,u±(P,r)|) dP 

Observe that we can pick the vector field H such that {H{P),N{P)) > 1 independent of P S dil. 
Hence applying the ellipticity condition (A3) to the left hand side of (471, the estimates (481 and 
(491 to the right hand side of (471 and also Cauchy’s inequality, we will finally obtain 

[ (|Vu±(P,r)p + |r|KP,r)|2)dP< / (|5,,u±(P,r)p + |u(P,r)|2) dP, 


which proves (431. 


The proof of (44) is a bit more involved and goes as follows. Given a vector field H and using 
0. an elementary calculation yields 


Y dk[{Hkdij - Hidkj - Hjdik) diudju] = ^ bjjdjudjU- 2t Yl^i {Hjdjuu), 


2=1 


where bij = dk{Hkaij — Hihkj — HjCLik)- Therefore, divergence theorem implies that 

n p 

(50) Y / nk{P)(Hk{P)'a,j{P)-H,{Pyakj{P)-H,{P)d^k{P))d,u^{P,T)d,u^{P,T)dP 
1.-1 JdQ. ^ ' 

P n n 

: / 2t ^Im (^Pi(X) i9ju(A,r) u(X,r)^ + ^ bij{X) diu{X,T) dju{X,T) 

i=l ij = l 


ijj,k—l ' 

= ± 


dx. 


Now since 


n p 

? {Hk{pya,yp) - H,{p)dkjiP) - Hypya^yp)) a,u±(p,T)a,u±(p,T)dP 

n p 

Y (- / nyp) \Hyp) dkyp) - nyp) yyp)] 5Y(p,r) 5 ,u±(p,t) dP 

^ Jon L J 


i^j^k—1 


i^j^k—1 


Uk (P) Hj (P) (P) - Hk (P) hy (P) (P, r) djU^ (P, r) dP 


ion 


[ nk{P) Hk(P) dij{P) diU^{P, t) djU^{P, r) dP) , 
JdQ ^ 
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equation ( |50| yields that 

n « 

XI / n-kiP) Hk{P)dij{P)diU^{P,T)djU^{P,T)dP 

n p 

= - V / nk{P)\H,{P)ak,{P)-Hk{P)'a,,{P)]d,u^{P,T)d,u^{P,T)dP 

Lj.fc=i 

n « 

- ^ / nfe(P) [i7,(P)a,fe(P)-i?fe(P)a,,(P)j a,u±(P,T)a,u±(P,T)dP 




/on 


P n n 

T / 2t ^Im ^Pi(X)i9ju(X, r)u(X, r)^ + ^ bij{X) diu{X,T) dju{X,T) 

i=i ij=i 


dX. 


Now we observe that for fixed j (resp. fixed i) the vector nfc(P)[Pi(P) dkj{P) — Hk{P) dij{P)] (resp. 
nk{P)[P[j{P) dik{P) — Hk{P) dij{P)] is orthogonal to the normal vector N. Therefore the first two 
terms on the right hand side of ( [5T| can be estimated by Jg^ |Vit^(P, r)| |Vt'u(P, r)| dP. Moreover, 
if we once again use the fact that {H{P),N{P)) > 1, then using the ellipticity condition (A3) and 
(511 we obtain 


(52) 


[ |Vu±(P,T)|2dP< [ |Vu±(P,T)||VTu(P,r)|dP 

JdQ Jdn 


Ida 

+ [ 2|r| |Vu(A,t)| |u(A,T)|dA+ / |Vi(A)| |Vu(A, t)P dA. 

dn± Jn± 


Now since (46) yields that 


[ \T\\VuiX,T)\\u{X,T)\dX< f |r|5 |Vu±(P,r)||u(P,r)|dP, 

Jq± Jon 


estimate (52) and Cauchy’s inequality imply that 


(53) [ |VM±(P,T)pdP< / (|VTu(P,T)|2 + |r||u(P,r)p) dP+ / |Vi(A)| |Vm(A, r)|2 dA. 

Jdn Jon ^ ^ Jn± 

As the reader can see, the estimate above contains an undesired term |VA(A)| |Vu(A,r)pdA 
which will be removed in what follows. To this end since A, verifies all the assumptions of Lemmas 
6.1 and |6.3[ using these lemmas we have (53) 


[ |Vu±(P,r)|2dP< [ (|VTu(P,r)p + |T|mP,T)|2)dP 
Jdn Jdn ^ ^ 

+ sup |Vu(A,r)|2dP+-^ / |Vu(A,T)pdA 

Jdnx(^~ts{P) " ’ Jn± 

< [ (|VTu(P,'r)|2 + |r||M(P,r)p)dP + (5”-i+“« [ sup |Vu(A,T)pdP 

Jdn^ ' Jdnxe-ysiP) 

where we have used the fact that \d^^u^\ < |Vu^|. Next picking £ sufficiently small, say e < £o(d) < 
1, we can absorb the term Jg^ |Vu^(P, r)pdP in the left hand side of the estimate above. 

This establishes (44). □ 


Now we are ready to state and prove our first Rellich-type estimates for the Fourier-transformed 
equation associated to the diffusion matrix A + A — A^‘^\ 


Theorem 6.5. Let A he a matrix satisfyi ng pr oper ties (Al) - (A4) 
with 0 < r < 1, are defined as is Lemmas 2.6 


and 


2.7 


and assume that A and A^'^f 
Set B := A -|- A — A^’’^ and assume that 


Cbu = 0 in LI. Then for all S S (0, do(''’)) G (0, £o(dA))> we have the following Rellich-type 

estimates 

(54) f (|Vu±(P,r)|2 + |r||u(P,r)|2^dP< / (|9,,u±(P,r)|2 + |u(P,r)p) dP, 
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(55) 


(56) 


I on 


[ \Vu^{P,T)fdP< [ (|VTw(P,T)|2 + |r||7z(P,r)|2 + ML^]dP 
JdQ Jon ^ e / 

+ <5"-i+“c f \{Wu)t{P,r)fdP, 

Jon 

\d,,U-{P,T)\^dP< [ +(p^^)|2+^n-l+»o|(V^)+(p,r)|2) dP. 

Jon ^ I 


where (Vu)J denotes the non-tangential maximalfunction ofVu. 


Proof. By definition, we have that B := A + A — A^'^'> and B satisfies properties (Al) - (A4) The 
idea is to transfer the Rellich inequalities (431 and (441, which are valid for the smooth matrix A, 
to the corresponding ones for B. 

Let mij be the matrix associated with A — A^'^\ and bij be that of B’s. Observe that by the 
construction of A^'^'> 

mij{X) = 0 when dist(X, 90) < ror, 
where rg is the diameter of O. Since Bbu = 0 in O, it is clear that 

n n 

C^u = —iru — — rriij) djU\ = di{mij dju) in O. 

i,j = l i,j=i 

Hence, we can decompose u = w + v; where = 0 in O, and 

Up n 

v{X,t):=-Y^ / dyT^{X,Y,T)mij{Y)dy^u{Y,T)dY =:'^Vij{X,T), 

_,40 ■ ■_ 1 


i,i=l ■ 


i,i=l 


where T^ denotes the fundamental solution T considered in Section]^ but here we emphasise its 
dependence on A. 

First we prove (54). We observe that applying (43) to w yields 
(57) 

[ (|Vu±(P,r)|2 + |r||u(P,r)|2)dP 
Jan ^ ' 

(|Vu;±(P,t)P + |r||u;(P,r)|2) dP + (|Vu±(P,r)|2 + |t| |u(P,r)p) dP 


< 


I an 


<f {\d,.w^{P,T)\^ + \w{P,T)\^)dP+ f (|Vu±(P,r)|2 + |r||u(P,r)|2)dP 

t/ 5r2 t/ dCl 

< [ (|9,,u±(P,r)p + |u(P,r)|2)dP+ f (\Xv^ {P,t)\^ + \t\\v{P,t)\^ + \v{P,t)\^) dP, 

where we have used that H = P on 90, and |9i/^u^| < |Vu^|. It remains to control the integrals 
involving the function v. 

Applying Lemma |3.8| (t), {ii) and (Hi) (taken with q = 1/2), we have that 

|Vu(F,r)| 


+ \'r\^^^\vij{P,T)\ + hj{P,T)\ < f 

Jn 


dY 


< 


nn{|V-P|>ror-} A ~ P\^' 

4 / |Vu(y,T)|dy, PS90. 
^ Jn 


Thus, 


(58) / (|Vz;±(P,r)p + |r|K,(P,r)|2 + |u.,(P,r)4dP<4^ / \VuiY,r)\^ dY. 

JdVt ^ ^ JQ, 

Now applying Lemma |6.l| to the term |VM(Y,T)pdy in the above estimate and inserting the 
resulting estimate in ( |57p , we arrive at 

[ (|Vu±(P,r)p + |r|KP,r)p)dP< (l+4:) / (|9.,u±(P, r)|2 + KP, r)^ rfP 
^00 ^00 


This proves (54). 
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We proceed with the proof of ( [55| . We observe that applying (441 tow yields 
(59) [ |VM±(P,r)pdP < [ |Vw±(P,T)|2dP+ [ \\'v^{P,T)\'^dP 

J dQ J dQ J dO, 


< 


Ian dP 

^gn-1+c.o f sup |Vw(X,r)pdP+ [ \\7v^iP,T)\^dP 
Jdnxei/.(P) Jan 


an 

<f (|VTu(P,T)|2 + |r||ir(P,T)p + M:^^)dP + ^”-i+““ / \{Vu)tiP,T)\‘^dP 
Jan ^ E / Jqq 

+ [ f\Wv^{P,T)\^ + \Wtv{P,t)\^ + \t\\v{P,t)\^ + 

Jan ^ 

^^n-l+ao f sup \Vv{X,T)\^dP. 

Jan xe-y/i(P) 


\viP,r)\^ 


Now using (581 and Lemma 6.1 we once again see that 


(60) 


[ (\Xv^{P,T)f + \XTviP,T)f + \T\\v{P,T)f + 
Jan ^ 


|p(P,r)p 


< 


< 


dP 

1 


[ (|Vi;±(P,r)|2 + |r|KP,r)p + Ml^)dp< 1 / \Xu{Y,r)\UY 

Jan ^ s ^ r Jq 


I an 
1 


I an 


(e'|Vw±(P,r)|2 + dP. 


Moreover using the support property of one has 


(61) 


\Xv,j{X,t)\ 


< 


|Vu(y,r)| 


dY. 


I an{ Y; dist(Y,an)>ror} |x-r|' 

If now P G dfl, X G jsiP) and dist(Y,9n) > rgr, then triangle inequality yields that 

|X - r| > ror - (1 + a),5 > C > 0, 


provided that we choose d < =■ do(r). Therefore choosing 6 < Sgir), estimate (611, Lemma 

and the Cauchy-Schwarz inequality yield 


6.1 


dP. 


(62) [ sup \Xv{X,T)\^dP< [ \Vu{Y,T)\‘^dY < [ (e'|Vu±(P, r)p + 

Janxe-fsiP) Jn Jan^ £ ^ 

Hence choosing e' small enough, say e' < £g(i5, r, e), we can absorb the term £'\Xw^[P,t)\^ dP 
on the right hand side of (601 and (621 in the left hand side of (59l. Thus, ( pOj ), (621 and (591 
together yield 

[ |Vu±(P,r)|2dP< [ (|VTw(P,T)|2 + |r||u(P,r)p + ^'^^)dP 

Jan Jan ^ e / 

+ <5"-i+“« [ \{Xu)t{P,T)\^dP 

Jan 


Finally, (56l follows by using the fact that \d^^u^\ < |Vu^| and applying first (55l and then 

(1^. □ 


6.2. Rellich estimates for the parabolic equation. Once again we start by some basic estimates 
for operators with smooth diffusion coefficients. However, for a Holder-continuous A we have the 
following lemma 


Lemma 6.6. Let A be a matrix satisfying properties (HI) 
and u{X, 0) = 0. Then for every £ > 0 


(H4), Assume that Cau = 0 infix (0, T) 



\Vu{X,t)\^ dX dt < [ [ (e|Vu±(P,t)p 

Jo Jan ^ 


HP,tr 

S 


dP dt. 
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Proof. We multiply equation Cau = 0 by u and integrate by parts. Thereafter we use the ellipticity 
assumption (A 3 ) and the assumption u{X, 0) = 0, which yield 


(63) 


± 


/o Jon 


u u* dP dt = 


> 


/ / 7 ciij diudjU + udtu ] dX dt 

Jo Jn± ) 

I I iXufdXdt+l I \u{X,T)fdX. 
Jo 2 Jq± 


Now, neglecting the last term in the right hand side of (631 and applying Cauchy’s inequality we 
obtain the estimate claimed in the lemma. □ 


Lemma 6.7. Let A S C°°{2Qq \ dft) be a matrix verifying (131, then for eaeh 5 > 0 one has 
\XA{X)\\Xu{X,t)\^dXdt<d'^-^+°‘° [ [ sup \Xu{X,t)\^dPdt 


10 


to Janxe-YsiP) 


+ 


1 


5(l-«o) 


'0 


\Xu{X,t)\^dX dt, 


where jsiP) was defined in Lemma 6.3 and is the parameter appearing in (13l. 
Proof. The proof is practically identical to that of Lemma [ 6 .3| hence omitted. 


□ 

To handle the Rellich estimates involving fractional derivatives of the solution of the parabolic 
equation, we would also need the following lemma 


Lemma 6 . 8 . Let Abe a matrix satisfying properties (Al) - (A4), Assume that Cau = 0 in fix (0, T) 
and M(7f, 0) = 0. Assume that H is a vector field on dft with {H,N) > 0. Then, for every 
e > 0, we have 

(64) [ [ {H,Xu{X,t))dtu{X,t)dXdt< [ [ (e\Xu^{P,t)\^ + dPdt. 

Jo Jn± Jo Jan ^ e ) 

Proof. Let v and w be solutions of Cau = 0 in fl x (0, T) and v{X, 0) = w{X, 0) = 0. Then, using 
the equation for v and w and the divergence theorem yield 

n p'J' p p'J' p 

(65) ^ / / Uij div{X,t) djw{X,t) dX dt + / / v{X,t) dtw{X,t) dX dt 

, .—I Jo Jfi^ Jo Jq^ 

= ± / f v{P,t)dvAW^{P,f)dP dt. 

Jo Jan 

Now if we take v = w = dI^^u (noticing that is also a solution of the equation) in (65) and 

use the ellipticity condition (A3), Lemma [2. 8 1 and the fact that n(A, 0) = w{X,{)) = 0 (to be able 
to eventually omit the second term on the left hand side), we obtain 

( 66 ) 


/7 

/o Jn^ 


\Dl'^Xu{X,t)\’^ dXdt< J dPdf 

- / (r \d..u^ip,t)\^dty^"(^ r \Di/\{p,t)\^dty^\p. 


Using (14) we observe that dtIi /4 = Furthermore, if we take v = and w = I 1 / 4 U we 

observe that these are still solutions of the equation and /i/ 4 u(X, 0) = 0. Therefore if we apply 
(651 to these solutions and use the ellipticity again we get 

r'^ r r r'^ 1/2 r'^ 1/2 

(67) / / \Df‘^u{X,f)\^dXdt< ( \d,,AU^{P,t)\‘^dt) ( / \Dl^‘^u{P,t)\‘^dt) dP. 

Jo Jn^ Jan ^Jo ^ ^ Jo ' 


A/4, 


tan '-Jo 

Thus ( |^ , (67) and Cauchy’s inequality yield 

fT 

( 68 ) 


to Jn^ 


(^\Dy\u{X,t)\^ + \D^/\{X,t)\ydXdt 
Jo Jan ^ £ ) 
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Finally using the fact that dtu = on the left hand side of (641, Lemma 

( 681 , we can conclude the proof of this lemma. 


2 .J 


and estimate 

□ 


Proposition 6.9. Let A be as in Lemmas 2.6 and assume that Cju = 0 m x (0, T) and u{X, 0) = 
0. Then 

[ [ \Xu^{P,t)\^dPdt< [ [ (|VTu(P,t)P + \Dy^u{P,t)\^ + dPdt 

Jo Jen Jo Jan ^ e J 

_^^n-i+ao f f gup \s/u{XA)\'^ dPdt, 

Jo JanxE-iaiP) 


for any 5 > 0 and small enough e € (0,£o(<^))- 
Proof. Following the proof of ( [44| we have 

rT 


E 

i,j,k—l ' 


0 JdQ 

pT 


= ± 


JO Jn± 


nk{P) (Hk{P)'a^j{P) - H,{P)~ajk{P) - Hj{P)~a,k{P)) d,u^{P,t) d,u^{P,t) dP dt 

r> n n 

2 ^ (^Hi{X)diu{X,t)dtu{XA)^ + ^ b^j{X) d^u{X,t) dju{X,t) 


i=l 


ij=i 


Hence, the same argument as in the proof of (441 yields 

rT r rT r 


JO Jan 


\Xu{P,t)\ dPdt 


< 


10 Jan 

.r 


\XTu{P,t)\\u{P,t)\dPdt- 


lo Jn± 


dX dt. 


(H, Xu{X, t)) dtu{X, f) dX dt 


/ / \XA{X)\\VuiX,t)\^dXdt. 

Jo Jn^ 

Now applying the Cauchy inequality to the first term on the right hand side of the above estimate, 
Lemma 6.7 and Lemma 6.6 to the second and Lemma 6.8 to the third, we obtain for every <5, £,e' > 0 

\Xu^{PA)\^dPdt< [ [ (\XTu{P,t)\^ + \u{PA)\^)dPdt 

Jo Jan ^ ' 


JO Jan 


_^^n-i+ao f f gup \Xu{X,t)fdP 
Jo Janxe-ysiP) 

rT 


dt 


(5(1 “o) 


(e|Vu=^(P,t)|^ 


\u{P,t)\^ 


dP dt 


JO Jan 


(e'\Xu^{P,t)\^ + dPdt 


Taking e and e' sufficiently small, say e < £o(<^) and e' < £q, we can absorb the terms containing 
|Vu^| in the left hand side of the estimate. This concludes the proof of the proposition. □ 

At this point we state and prove our second Rellich-type estimates for the parabolic operator 
associated the the diffusion matrix B. 

Theorem 6.10. Let B be as in the statement o/Theorem |6.5| and assume that Cbu = 0 m fix (0, T) 
and u{X,0) = 0. Then for all sufficiently small 6 € (0, <5o(7’)) and £ G (0,£o((5)), we have 

[ [ \Xu^{P,t)\^dPdt< [ [ (|VTu(P,f)p + |Dy\(P,t )|2 + MiV)P\^p^^ 

Jo Jan Jo Jan ^ £ J 


rn—l+ao 


\{Xu)t{PA)\^ dPdt. 


Jo Jan 


Proof. We are going to follow the scheme of the proof of Theorem |6.5[ So unless otherwise stated, 
we use the same notation as there. We decompose u = w + v, where this time C^w = 0 in 17 x (0, T), 
and 

n „j’ „ n 

■u(A,t) = -Y / dy.T^{X,Y,t-s)m,j{Y)dy.u{Y,s)dY ds v,j{X,t). 

.j^T^Jo Jn 


ui=i 
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First we apply Proposition |6.9| to the function w to get 

pT 

(69) 


10 Jdn 


[ \Vu^{P,t)\‘^dPdt< f f \Vw^{P,t)\‘^dPdt+ f f \Vv^{P,t)\‘^dPdt 
Jon Jo Jon Jo Jon 

[ f (\VTw{P,t)\^ + \Dl/^w{P,t)\^ + \wiP,t)\AdPdt 
+ ^n-i+ao f f sup \\/w(x,t)\^dPdt+ f [ \\7v^{P,t)\^dPdt 

Jo JdO.X^jsiP) Jo J dQ 

,.'T 

(|VTM(P,i)P + \Dl^^u{P,t)\^ + |u(p,i)p) dPdt 


Jo Jan 

+ 

V 7 

Jo Jan 


lo Jan 


\{Wu)4X,t)\^ dPdt 


(|Vz)±(P,i)p + \Dl/^v{P,t)\^ + HP,t)\^) dPdt 
^gn-l+ao j f sup \\7v{X,t)\^dPdt. 

Jo JdQ XGy^(P) 


>0 JdQ.X^jsiP) 

Let i, j = 1, ... ,n. Now Lemmas 3.3 3.4 (i); and 3.5 (ii) yield 

|Vn±(P,t)| + \Dy\,{P,t)\ + \v.,{P,t)\ < r [ 

Jo Jn 


< 


0 v/nn{|Y-P|>ror} 
1 


|v»(Ls)l 

\Y-P\n+2 


dY ds 


[ [ \Xu{Y,s)\dYds, {P,t)eST- 

’’’ Jo Jn 


Hence nsing Lemma |6.6| we obtain 

rT 


10 Jan 


(|Vn±(P,t)|2 + \Dy\{P,t)\^ + |n(P,t)p) dPdt 


< 


r^^+'^Jo Jn 


\^u{Y,s)fdYds<^[ [ (e|Vu±(P,t)|2 + 77V7)dPdt. 

^ Jo Jan ^ £ / 


Moreover, the same reasoning as in the proof of Theorem 6.5 and Lemma 6.6 yields for <5 < (5o(^) 
that 


f f sup \Xv{X,t)\'^ dP dt < ( f \Xu{Y,t)\‘^ dY 
Jo Janxeis(P) Jo Jn 

cT r , 

(s\Xu^{P,t)\ 


HP,tr 


dP dt. 


Jo Jan 

To conclude we only need to pick e small enough such that we can absorb the terms corresponding 
to ejVu^p in the left hand side of (691, which finally yields 

[ [ \Xu^{P,t)\^dPdt< [ [ (\XTu{P,t)\^ + \Dl^\{P,t)\^ + \u{P,t)\^)dPdt 

Jo Jan Jo Jan ^ J 

_^^n-l+ao f j \{yu)t{X,t)\'^ dPdt. 

Jo Jan 


□ 


7. Invertibility of operators associated to the layer potentials 


7.1. Invertibility of BSI. The main result of this Section is the proof of the invertibility of BSI. 
To achieve this we shall start by showing the invertibility of the Fourier-transformed layer potentials 
associated to suitable pieces of the diffusion matrix A. 


Lemma 7.1. Let B he the diffusion matrix in Theorem 6.5 Then for t 0 the operators ±1/2 ± 
defined in Remark 5.1 verify the estimate 


ll9lU.(»«. S (l + l/i) II ( ± 5 + 
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Proof. We only give the proof for 1/2 + (K^)*, since the proof for —1/2 + (K^)* is almost identical. 
Let g G L^{dil) and set u{X, t) := S'^(g)(X). By the jump relation in Proposition 5.9 we can write 

(70) 3=Q + - ( - ^ + - d^,u-. 

From now on we drop the subscript B in the normal derivatives. 

Next, since u{X,t) = S'^{g){X), we have that Lbu = 0 in fl and for 6 < 6Q{r) (see the proof of 
the second part of Theorem | 6 . 5 [ ), Rellich estimate (56 1 and estimate (38 1 (i.e. the boundedness of 
(VS(^)J) altogether imply that 

(71) [ \d,u-{P,T)fdP< [ (\u{P,T)f + \d,u+{P,T)f + S^-^+‘^°\g{P)\^) dP. 

Jdn Jdn ^ ' 

Now Lemma 6.1 yields (observe that t ^ 0) 


/n± 


^[Y.T)\^dY< f (^|5,^u±(P,r)|2 + 

Jdfl ^ +1 


2 , HPxW 


and 


[ \Xu{X,T))\^dY< [ (e\d,,u^{P,T)\^ + 

Jn± Jdn ^ 

Therefore adding these two estimates we obtain 


e T 




dP, 


dP. 


(72) (|Vu(X,r)|2 + |u(X,r)|2) dX < (e(l + ^)|a,,u±(P,r)|2 + ^(l+^)|u(P,r)p) dP. 


On the other hand (721 and Lemma 6.2 yield 


[ \u{P,T)\^dP<C [ (\Xu{X,T)\^ + \u{X,T)\^)dX 

Jdn Jn± ^ ^ 

<Ce(l + ^) f \d,u+{P,T)\^dP+-{l + ^) [ \u{P,T)\^dP. 
^ Jdn e \t\ Jdn 


Thus taking e = 20(1 + l/|r|)) we obtain 


(73) 


[ \u{P,T)\^dP<(l + ^Y [ \d,U+{P,T)\^dP. 
Jdn ^ PK Jdn 


Thus combining (0 and ( [T^ we have 

(74) f \d,U+iP,T)\^dP+ f \d,U-iP,T)\^dP<(l+^Y f \d,U+iP,T)\^dP 
Jdn Jdn ^ PK Jdn 


+ 5”-i+““ f \g{P)\^dP. 
Jdn 


Therefore, using ( [70| and ( |74| , we obtain for S < So{r) 
Il5llz,2(an) < \\d^u'^\\L2(dn) + \\diyU~\\L2{dn) 


< 


4 + y)ll(^<»^-)> 


L 2 (an) 


+ C6 ^~^\\g\\L^(dn)- 


Now if we take <5 < min{(5o('r), (^)^7(" i+“o)| then 

IlsIU+an) ^ ||9i.u+||i2(9Q) + \\d^u~\\L2(dn) - (l + j/rj) (2 + (^b)*)^ 
which is the desired estimate. 


L 2 (on) 


□ 


The final step here is to show the invertibility of the Fourier-transformed BSI associated to the 
matrix A. 

Lemma 7.2. Let A satisfy (Tl)| - |(Pl4)| defined as in Lemma 2.7 Then for r 7 ^ 0 the operators 
±1/2 + (K(^)*, are invertible in Lf{dLl) and satisfy the estimate 

ll9P=,», £ (1 + /j) II (± ^ 
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Proof. As before, we write A as B + C where B := A + A — and C := — A. Then we 

observe that 


+ {^bT + 


and Lemma [TT] yields that 

ll-?llL^(ao)< (l + ^)||(±^ + (K^,r)5 


L^{dn) 


gGL^idil). 


Moreover since ±1 + (K^)* = ±| + (K^)* — (K^)* + we have 


L 2 (an) 


L^{dn) 


\\9\\LHon)<{l + ^^\\{±l + {^Br)9 

< (l + ^) I ( ± i + (K^)* - (K^)* + (K^)*)g 

< (l + ^) I ( ± i + (K^)*)g|U.(an) + (l + ^) || H^a)* - (K^)*)5 

Now using Theorem 5.7 (estimate @) and Lemma |2.7[ we have 

((K^r - (K^D-? < \\A-B\\l/^gU.^ea) < \\C\\l/^g\\man) 

L^(ar 2 ) 

< \\A^^^ - AWU^gWman) < cr^^^^gWman), 


L^{an) 


where the last estimate follows from Lemma 2.7 Therefore choosing r so small that c(l+l/|r|)r“°/^ < 
1 /2, we will have 

iMi..,»,<(i + /|)||(±i + Kr)4.,^„. 

This shows that the operator 1/2 + (K^)* : Lf{d^) — LfidPi) is injective and has closed range in 
L^{dVi). To conclude, we need to show that the range of 1/2 + (K^)* is equal to Lf{dVl). We write 

(75) 


(2 +K^) - (2 + (n - (Kr) + {K-^a)‘ 


First we observe that Lemma 3.8 (iv) and Lemma 2.3 yield that is a compact operator. 


hence (K )4 — K^)* is also compact. Now we claim that the operator — (K));)* is compact as well. 
Indeed, let TAiX,Y) be the fundamental solution of the elliptic operator V • (AV). Then as was 
shown in m Theorem 2.6] one has 


where 


Eo{X,Y) = 


Ta{X, Y) = Eo{X, Y) + EiiX, Y), 
Cn 


(det A(F))i/2(A-i(y)(X -Y),X- F)("-2)/2 
and |Vjcifi(A, y)| |A — y|i-"+“/2^ where a is the Holder exponent of A. 

Now since the integral kernel of the operator (K)];)* is the normal derivative of r^(A, Y), estimates 
above yield that the integral kernel of K° — (K]]^)* is of the form 

(Vv((A-i(A)(y -X),Y- X)),A-\X)N{Y)) {Xx{{A-\Y)iX -Y),X- Y)),A-AY)N{X)) 
(det A(A))i/2(^-i(X)(y -X),Y- X)"/2 (det A(r))i/2(A-i(y)(x -Y),X- r)V2 ’ 

plus a kernel, which in light of the estimate for XxEi{X,Y) above, gives rise to a compact operator. 
This implies that the integral kernel of K/j — (K/j)* can be estimated by |A — Tmodulo a kernel 


which gives rise to a compact operator. Hence, using Lemma 2.3 we conclude that IK° — (IK°)* is 
compact. Thus, the operator on the left hand side of ( [7^ is compact and since 1/2 + (K^)* is 
injective. Lemma 2.5 yields the surjectivity of 1/2 + (K^)*. 

□ 

Now we are ready to state and prove one of main results of this paper, which concerns the invertibility 
of the boundary singular integral operator which is the backbone of the solvability results in Section 


Theorem 7.3. Assume that A verifies (Al) - (A4) Then, the operator ±1/2 ± defined in 
Remark 4.2, is invertible in L'^(St). 
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Proof. We only show the invertibility of 1/2 + since the invertibility of —1/2 + JC*^ is done in 
a similar way. First we show the injectivity of 1/2 + IC\. To this end suppose / S LP'{St) is such 


4.10 


that (1/2 + /C/i)/ = 0. Then if u{X,t) = Sf{X,t) then the jump relation in Proposition 
yields that d^u~^ = 0. Now using ([6^, it follows that u = 0 in x (0,r) and hence u = 0 in St- 


Furthermore, (63) also yields that u = 0 in x (0,T), and therefore d^u~ = 0. Thus once again 
the jump relation in Proposition |4.10 yields that / = 0 and hence (1/2 + /C/t) is injective. In order 
to prove the surjectivity of (l/2 + /C(/) we follow the strategy of Shen in |141 Theorem 4.1.18]. Take 
an arbitrary / € L^{St) and extend it to 917 x K by setting 


if 

hiP,t):= {-f{P,t-T) if 
0 if 


0 < t < T, 

T <t<2T, 

2T < t or t < 0. 


This function is in L^(917 x M) n L^(917 x M), and therefore its Fourier transform in t, i.e. h{P,T) 
is continuous in r. Furthermore h{-,T) G L^(dfl) for each r G K and h{P,0) = 0. Using Plancherel’s 
theorem it follows that ||/i||i2(aQxR) = C'l|/»llL2(aoxK) ^ WIUlhSt)- Now if r ^ 0 then using the 
invertibility of 1/2 + (K^)* on L^(917) established in Lemma 7.2 there exists ip{P,T) G Lf{dVl) 
such that 

(76) (l/2 + (K^)*)<p(P,r) = MP,r). 

For T = 0, thanks to the fact that h{P, 0) = 0, we can set (f{P, 0) = 0 and hence produce a solution 
to (761 for all values of r. We will show that 


(77) ll7^llL2(9nxR) < C't||/||l2(St)- 
Indeed, Lemma [7.1] yields that 

(78) IIV5(GT)llL2(an) < 


1 


(l + ■r:T)ll^(G'^)llL2(an). 


< 


other hand for |r| < 1, 


Therefore J^^^^Jg^\<f{P,T)\‘^dPdT < \\h\\l,^g^^ 

using the fact that h(P, 0) = 0 and applying the mean value theorem and the Cauchy-Schwarz 
inequality we deduce that 


\h{P,T)\<\T\ t|/i(P,7)|dt < Cr|T|{ / |/(P,t)pd7}^/^ 

^0 Jo 

This together with ( [7^ yield that Jqq | 7 ’(U, t)P dP dr < Ct Il/lli 2 (s,j,)- This establishes ( [77| . 
Now let ip{P,t) G L^{dfl X K) be such that its Fourier transform in t is equal to (p{P,T). The for 
almost all (P, t) G 917 x M we have 

l:i’{P,t) + [ [ X\{P,Q,t - s)'ip{Q,s) dQds = h{P,t). 

2 J_aa Jon 

Furthermore, since h{P, t) = 0 for t < 0, an argument similar to that of the injectivity of 1/2 + 
shows that 'tjj{P,t) — 0 when t < 0. Now, since for (P, 7) G St, h{P,t) = f{P,t), we have for almost 
every (P, t) G St that 

i'0(P,7)+ / [ X\{P,Q,t-s)ip(Q,s)dQds = {1/2 + X*)ip{P,t) = f{P,t). 

2 Jo Jdn 

Therefore 1/2 + is also surjective and hence because of its injectivity, invertible. □ 


Corollary 7.4. Assume that A satisfies the conditions (Tl) - (44) Then, the operators ±1/2+ /C^ 
are invertible in L‘^{St). 

7.2. Invertibility of the SLP. For the solvability of the regularity problem we would also need 
the following invertibility result. 

Theorem 7.5. Assume that A satisfies the conditions (41)| - (44), Then, the single-layer operator 
Sa is invertible from L'^{St) to H^’^^'^{St). 

Proof. First we show that Sa is injective. To this end take / G L'^{St) and set u := Ssf, with b 
as in Theorem 6.10 Then Cbu = 0 and it(X, 0) = 0. Moreover the jump relation in Proposition 
|4.10| yields 
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f={l + (^b)*)/ - ( - ^ + (/Cb)*)/ = 


- dvr,u 


Therefore the Rellich estimates in Theorem 6.10 (for 5 < ^o(?')) and the boundedness result (321, 
yield that 


II/IIl2(St) ^ W^vb^'^Wl'^^St) + 11^ “ IIl2(St) 

< |1Vtm||l2(St,) + ||T>y^M||L2(s^) + 111x11^2(3^) + 5 =^|i(V5B/)*||L2(s^) 

< ||u||bi. 1/2(3^) +(5 5 ||/||L2(Sr)- 

So we arrive at the invertibility estimate 

II/IIl2(St) ^ II‘5b/||bi.i/2(3j,) +(5 5 ||/||b2(3^). 

However, since Sb = Sb — Sa, we also have 


(79) ||/||l2(3,j,) < ||5a/||bi,i/2(3,j,) + \\{Sb - SA)f\\m,i/2(^ST) +^ II/I|l2(St) 

= II‘5a/|1bi,i/2(3j,) + ||V(5b - 5yl)/||i2(3,j,) + \\Dy^{SB “ 5yl)/|| 32 ( 3 ,^,) 

+ II(5b - 5a)/||l2(3.j,) + 6 ^||/||l2(St). 


Now using Theorem |4.14| and Lemma [2.7| we can show that the last three terms on the right hand 
side of (79l are bounded by ||H — B\\ii‘^\\f\\B2(^ST) ^ ^““^^II/IIl 2 ( 3 ,j,). So taking r small enough we 
can absorb the second, third and fourth terms on the right hand side of (791 into the left hand side, 
and thereafter choosing 5 in a suitable manner, we can do the same with the fifth term on the right 
hand side of (791 and finally obtain 


(80) II/IIl2(3,) < 115^/11 

Thus if IX = 0 then ( |80| yields that / = 0, and therefore 5^4 is injective and its range is closed 
in (5 t). The surjectivity of the single layer potential is quite standard and is handled in 

many papers, see e.g. (2 Theorem 5.6 ] or m Theorem 4.2.1]. We would like to mention that, the 
main part of the surjectivity-proof, hinges on two main steps; one is an approximation argument 
and the existence of regular solution to the Dirichlet problem for equations with Holder-continuous 
coefficients on smooth domains, see [5] . The other is the existence and uniqueness of the solution 
to the Neumann problem, which will be established in Theorem |8.2[ 

□ 


8. Solvability of initial boundary value problems 

Given all the information and estimates which have been gathered in the previous sections, one has 
the following three theorems which state the well-posedness of the DNR problems. The proofs of 
these statements are completely standard and can be done in the same way as e.g. [S], hence we 
omit the details. 


8.1. Solvability of the Dirichlet problem. 

Theorem 8.1. Let T > 0, fl be a Lipsehitz domain and A a matrix verifying properties (Tl)| - [(T4)[ 
If f & L‘^{St), then the function u := I?(—1/2 -|- is the unique solution of the following 

initial Dirichlet boundary value problem 

dtu{X,t)-^x ■ {A{X)Vxu{X,t)) =0 in Ll x {0,T), 


x(X, 0 )= 0 , a.e.XGD, 


and 


lim u{X,t) = f{P,t)i {PA) ^ St, 

x^p 

XG7+(B) 


k?llL2(3,r) < ||/||l2(3,j,). 


Proof. The proof goes along the same lines as [2 Theorem 2.3]. 


□ 
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8.2. Solvability of the Neumann problem. 


Theorem 8.2. Let T > 0, LI be a Lipschitz domain and A a matrix verifying properties (Tl) - 

~^f is the unique solution 


(A4), If f G L^{St) with f = 0, then the function u \= Sa{^/‘ 2- + 
of the following initial Neumann boundary value problem 


dtu{X,t)-'L'x ■ {A{X)\/xu{X,t)) =0 in nx{0,T), 


uiX,0)=0, a.e.XGLl, 


di,^u{P,t) = a.e. {P,t) & St, 


and 


Wu^Wl^St) + ll(Vu)f + \\{Dy^u)t\\L2(ST) ^ II/IIl^(St)- 

Proof. The proof goes along the same lines as [5] Theorem 2.4]. 

8.3. Solvability of the regularity problem. 


□ 


Theo rem 8.3. Let T > 0, LI be a Lipschitz domain and A a matrix verifying properties |(A1)| - 
(A4 ). If f G then the function u := )~^ f the unique solution of the following 

initial regularity boundary value problem 

dtu{X,f)-Vx ■ {A{X)Vxu{X,t)) =Q in Ll x {0,T), 


u{X,0)=0, a.e.XGLl, 


lim u{X,t) = f{P,t), a.e. {P,t) G St, 

JC — y}-^ 

xej+{P) 


and 


Proof. The proof goes along the same lines as Theorem 3.4]. 


□ 
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